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ABSTRACT 

In this paper we present high resolution VLT/FLAMES observations of red giant stars in the 
massive intermediate-age Large Magellanic Cloud star cluster NGC 1846, which, on the basis of its 
extended main-sequence turn-off (EMSTO), possesses an internal age spread of ~ 300 Myr. We 
describe in detail our target selection and data reduction procedures, and construct a sample of 21 
stars possessing radial velocities indicating their membership of NGC 1846 at high confidence. We 
consider high-resolution spectra of the planetary nebula Mo-17, and conclude that this object is also 
a member of the cluster. Our measured radial velocities allow us to conduct a detailed investigation 
of the internal kinematics of NGC 1846, the first time this has been done for an EMSTO system. 
The key result of this work is that the cluster exhibits a significant degree of systemic rotation, of 
a magnitude comparable to the mean velocity dispersion. Using an extensive suite of Monte Carlo 
models we demonstrate that, despite our relatively small sample size and the substantial fraction of 
unresolved binary stars in the cluster, the rotation signal we detect is very likely to be genuine. Our 
observations are in qualitative agreement with the predictions of simulations modeling the formation 
of multiple populations of stars in globular clusters, where a dynamically cold, rapidly rotating second 
generation is a common feature. NGC 1846 is less than one relaxation time old, so any dynamical 
signatures encoded during its formation ought to remain present. 

Subject headings: globular clusters: individual: NGC 1846 — Magellanic Clouds — stars: kinematics 
and dynamics 



1. INTRODUCTION 

1.1. Multiple populations in Galactic globular clusters 

One of the outstanding problems in modern astro- 
physics concerns the formation of globular clusters. Long 
thought to constitute a simple, homogeneous class of 
object, each one consisting of stars of a uniform age 
and elemental composition, it is now recognised that 
these systems harbor multiple stellar populations dis- 
playing a wide varie ty of unexpected characteristics (see 
iGratton et alll20Tl for a review . 

All Galactic globular clusters for which large samples 
of members have been studied spcctroscopically at suf- 
ficiently high resolution and signal-to-noise are found to 
be comprised of stars exhibiting a characteristic chemical 
signature - a strong anti-correlation between the abun- 
dances of the light elements O-Na, as well as C-N and 
Mg-Al in many cases, e ven while remaining ho mogeneous 
in iron content (see e.g. JCarretta et a l. 2009a, hQ). Since 
this pattern is observed in stars on the main sequence as 
well as in red giants, it must be primordial rather than 
being the result of nucleos ynthesis and mixing within the 
observed giant stars (e.g. JGratton et al~l l2004) . The pat- 
tern is also seen in ancient globular clusters in nearby 
galaxies - for example, in the Large Magellanic Cloud 
dHill et al.l [20001: I Johnson et all 120061: IMucciarelli et al.l 
2009) an d the Fornax dLetarte et al.ll2006D " and Sagittar- 
ius (e.g., ICarretta et al.l I2010allbl ) dwarfs - but not ob- 

• Based on observations obtained at the European Southern 
Observatory Very Large Telescope, Paranal, Chile, under pro- 
gramme 082.D-0387. 



served in ol d open clusters (wi th the recent exception of 
NGC 6791. lGeisler et al.ll2012h or for the vast majority 
of the Milky Way field h alo (see e.g.. iMartell fc Grebell 
120101 : IMartell et al.ll201ll and references therein) , imply- 
ing that it is a product specifically linked to globular 
cluster formation processes. 

Beyond the anti-correlated light element abundance 
variations many globular clusters exhibit various addi- 
tional levels of inhomogeneity, such as split main se- 
quences or multiple sub-giant branc hes on their color- 
magnitude diagrams (CMDs) (e.g . . iBedin et al.l l2004t 
Piotto et all 120071: iVillanova et all [2007; M ilone et all 

2008h . internal dispersions in the abundan ce of 

iron or other elemen ts (e.g.. |Yong fc Grundah] 120081: 
| Da Costa et al.l 120091: [Ferraro et al I l2009t iCohen et al.1 
120 10t IMucciarelli et al.ll2012f) . and /or subpopulatiqns en- 
hanced in helium ( e.g.. iNorrisI [200l iPiotto et al.l 120051: 
iDupree et al.ll201ll iPasauini et al.ll2011l ). 

The overall picture is of a huge wealth of complexity 
that poses serious challenges for models of globular clus- 
ter formation and evolution. The pervasive light element 
anti-correlations require material processed at high tem- 
peratures via proton capture reactions. At T > 2 x 10 7 
K the CNO and NaNe cycles serve to alter the car- 
bon, nitrogen, oxygen and sodium abundances, while 
at somewhat higher temperatures the MgAl cycle also 
becomes active, leading to enhanced aluminium and re- 
duced magnesium abundances. In the presently favored 
model, material processed in this way in a first gener- 
ation of stars pollutes or forms a central accumulation 
of gas in a young globular cluster, and a second gen- 
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eration of stars bearing the characteristic light element 
sign ature is subsequently formed from this reservo ir (see 
e.g.. lD'Ercole et alJl200a IConrov fc Spergelll2011fl . Note 
that because the iron content in most clusters is ob- 
served to be homogeneous, the gas should typically not 
have undergone supernova enrichment. Leading sugges- 
tions for the sites of the high temperature processing are 
inter mediate-mass asymp totic giant branch (AGB) stars 
(e.g.. lVentura et aT1l200ll) and fast - rotating massive stars 
(FRMS: e.g.. lDecressin et al.ll2007[ ) . In both cases helium 
enhancement may also be introduced to the gas reservoir 
as a result of main sequence hydrogen burning in the 
polluting stars. 

A generic prediction of this model is that the second 
generation of stars in a globular cluster ought to be more 
centrally concentrated than the first generat ion, and this 
seem s to have been observed in some cases (jLardo et al.1 
1 2 1 lh . However, the model also suffers from a number of 
difficulties, chief among which is that the second genera- 
tion stars in globular clusters require significantly more 
gas for formation than can have been available based on 
the presently-observed numbers of first generation mem- 
bers. Ideas invoked to circumvent this issue include the 
accretion o f large amounts of prist ine interstellar ma- 
terial (e.g.. IConrov fc Spergell |20TTI ). a first generation 
~ 10—100 times more massive than p r esent-day gl ob- 
ular clusters (e.g., iD'Ercole et all [20081: iBekkil [201lh . a 
first generation with a top-heavy initial mass function 
(e.g.. lD'Antona fc Caloill2004t IBekki fe Norrisll2006[ h or 
a variation in which globular clusters are formed as part 
of initially larger systems such as low- mass dwarfs (e.g., 
IBekki et al.ll2007|) . Beyond this difficulty, it is also un- 
clear how clusters with internal dispersions in iron-peak 
elements, or that include populations with very high he- 
lium abundances (Y ps 0.4) fit into the model. 

The unavoidable conclusion from the scenario outlined 
above is that the formation of individual globular clusters 
must have spanned a period of tens to hundreds of Myr, 
depending on the nature of the stars responsible for the 
high temperature processing. With presently available 
facilities, however, we are unable to directly resolve age 
differences of this magnitude given that Galactic globular 
clusters are typically ~ 12 Gyr old. In addition, since 
the time-scale for dynamical relaxation in the majority 
of Galactic glob ular clusters is much shorter than the 
cluster age (see iHarrisi 1 1996T ) any detailed information 
imprinted on the internal kinematics of these systems as 
a result of the formation process will have long since been 
seriously diluted or possibly erased altogether. 

1.2. The role of peculiar Magellanic Cloud clusters 

Star clusters in the Large and Small Magellanic Clouds 
(LMC and SMC) offer an important new piece of this 
puzzle. These two galaxies possess extensive systems of 
clusters spanning the full age range ~ 10 6 — 10 10 years; 
many of the members of these systems are comparable 
in mass to present-day Galactic globular clusters lying 
at o r below the peak of the luminosity function (see 
e.g.. iMackev fc Gilmorel l2003alTbh . Using images taken 
with the Advanced Camera for Surveys (ACS) onboard 
the Hubble Space Telescope (HST), we recently demon- 
strated that several rich intermediate-age clusters (r ~ 
1.5 — 2 Gyr) in the LMC di s play markedly unusual CMD s 
(|Mackev fc Brobv Nielse"n1 120071 IMackev et all I2008af> . 



More specifically, we found that while the CMDs for 
NGC 1783, 1806 and 1846 have very narrow main se- 
quences and red giant branches, the main-sequence turn- 
offs for these clusters are much broader than can be ex- 
plained by the photometric uncertainties (see Figure [1} . 
After considering various possibilities for the origin of 
these extended main-sequence turn-offs (EMSTOs) , such 
as confusion due to unresolved binary stars or field con- 
tamination, we concluded that the simplest viable inter- 
pretation of our observations is that each of the three 
clusters is comprised of two or more stellar populations 
with very similar iron abundance but spanning an age 
interval of « 300 Myr. 

Subsequent work has reinforced this interpretation and 
revealed that an EMSTO is apparently not an unusual 
feature for intermediate-age Magellanic C loud clusters. 
Photometr ic analysis of HS T ima ging b y iMilone et al.1 
(|2009fl and lGoudfrooii et all (|2009L l2011al ) demonstrated 
that, of the 16 intermediate-age LMC clus ters with suit- 
able data, 11 posses s EMSTOs ( see a lso IMackev et al.l 
12009ft . In addition, iGlatt et al.l (12008ft discovered that 
the intermediate-age SMC cluster NGC 419 possesses an 
EMSTO. Each of these studies concluded that an inter- 
nal age spread was the most likely explanation for these 
features, with the full sample of EMSTO clusters encom- 
passing the range ~ 150 — 500 Myr for this spread. 

Additional s uppor t has come from several directions. 
IGirardi et al.1 (|2009ft noted the presence of a dual red 
clump in NGC 419 an d in a number of the rich LMC 
EMSTO clusters, while iRubele et all poTol 12011ft used 
the complete observed CMDs for NGC 419 and NGC 
1751 to reconstruct their star-formation histories - find- 
ing that these may have spanned an incredible ~ 700 Myr 
(out of a mean age of ~ 1.5 Gyr) in NGC 419 , and ~ 460 
Myr in NGC 1751. IRubele et all IpOloL l20ll noted that 
in their best-fitting models the dual red clumps in both 
NGC 419 and NGC 1751 arise as a direct result of t he age 
dispersions in these systems. iBastian fc de Minkl (2009) 
posited that rather than reflecting an internal age spread, 
a cluster EMSTO might instead arise if a wide range in 
ste llar rotation is presen t at the MSTO. However, models 
bv IGirardi et al.l ()2011f ) demonstrated that the effects of 
stellar rotation fail to reproduce th e EMSTO morphol- 
ogy; in addition, as emphasised by IRubele et al.l (|2010t 
1 2 llf) , this scenario does not seem able to account for the 
dual red clumps observed in NGC 419, NGC 1751 , and 
some other EMSTO clusters. IGirardi et all (|2011l ) also 
tested the effects of star-to-star variations in the degree 
of convective overshooting on a cluster CMD, but found 
that this too failed to accurately reproduce the charac- 
teristic EMSTO shape. 

The global properties of EMSTO clusters provide 
some clues as t o how these systems may arise. 
IConrov fc Spergel (|2011h made the simple observation 
that it is only the intermediate-age clusters with present- 
day ma sses greater t han w 10 4 M Q that exhibit EM- 
STOs. Keller et al.l (|201lD further observed that the 
known EMSTO systems are, without exception, the most 
diffuse, spati ally extended clust ers f or their age. In the 
framew ork of lElson et al.1 (|1987l ) and Ma ckey fc Gilmorel 
(|2003airbT ). Magellanic Cloud clusters exhibit an increas- 
ing spread in size (defined by e i ther t he core or half-light 
radius) with age; iKeller et al.l (|2011f) showed that at in- 
termediate ages the EMSTO clusters all fall towards the 
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upper envelope of this distribution, while clusters with- 
out an EMSTO app ear systematically more compact. 
iMackev et all (|2008b[ ) used detailed N-body models to 
explore the origin of the increasing spread in Magellanic 
Cloud cluster sizes with age. They found that if very 
young massive clusters are formed as compact dense sys- 
tems, as is seen to be the case in both the LMC and 
SMC, the only viable way for them to evolve along the 
upper envelope of the observed age-size distribution on 
a < 1 Gyr time-scale is due to mass-loss from stellar 
evolution if they were initially highly mass segregated - 
that is, with the highest mass stars preferentially located 
towards the cluster centerfl This in turn implies that, 
in addition to being the most massive intermediate- age 
clusters in the Magellanic Clouds, EMSTO systems were 
also probably the most strongly mass segregated clusters 
at early times. There is clearly a link between cluster 
mass and structure, and the presence of an EMSTO. 

There are two additional observations of relevance. If 
a cut is made on the CMD across a cluster's EMSTO in 
a direction perpendicular to the locus of the upper main 
sequence, then groups of "younger" and "older" MSTO 
stars may be defined. These groups have di stinct proper- 
ties. First, as noted bv lMilone et al.1 fj2009f ) . the younger 
population consists of at least as many stars as the older 
population - in fact the ratio is typically more like 2:1 
in the richest clusters. Sec ond, as described in detail 
bv IGoudfrooii et all (|2011bt ). the younger population is 
often more centrally concentrated than the older popu- 
lation. The strongest difference between the two con- 
centrations is seen in the clusters possessing the largest 
estimated initial escape velocities (i.e., v esc > 15 kms -1 
at an age of ~ 10 Myr); systems in which v csc < 10 
kms^ 1 do not show this difference in central concentra- 
tions, or are not seen to possess an EMSTO at all (note 
that the limiting velocity here is comparable to the veloc- 
ities of winds from FRMS or AGB stars). This picture is 
fully consistent with the idea discussed above that EM- 
STO clusters were both the most massive and the most 
strongly mass segregated clusters at the time when they 
were formed. 

As speculated by a number of authors 
(IConrov fc Spergell 12011 IGoudfrooii et~aT1 I2011bt 
iKeller et all 1201 ID . the properties of EMSTO systems 
suggest a formation process remarkably similar to that 
inferred for the multi-population Galactic globular 
clusters - specifically, that prolonged star formation has 
occurred at the centres of these objects because their 
masses and initial structures allowed the retention or 
accumulation of a suitable reservoir of gas at the bottom 
of the cluster potential well. The observed internal age 
spreads of several hundred Myr in EMSTO clusters, and 
the apparently minimal dispersions in iron abundance 
inferred from their narrow RGB sequences are both con- 
sistent with the scenario invoked for Galactic globular 

2 Mackcv ct al. ( 2008b) found that a retained population of stel- 
lar mass black holes can also cause a cluster to move towards the 
upper envelope of the age-size dist ribution, bu t on l onger time- 
scales than 1 — 2 Gyr. Note also that Elson et al. ( 1989) found that 
a flat, or top-heavy IMF can lead to cluster expansion along the 
upper envelope of the age-size distribution; however such clusters 
rapidly become unbound (after just a fewxlO 7 yr). In addition, 
young Mag ellanic Cloud clusters are seen to have quite normal 
IMFs (e.g.. lKroupall200ll : lde Griis et al.ll200l ). 



clusters in which much of this reservoir comes from the 
slow winds of a first generation of AGB stars. Further- 
more, the number ratio of younger to older groups of 
stars in EMSTO clusters is comparable to that seen for 
the two generations in Galactic globular clusters. Note 
that it has not been clearly assessed how early ideas for 
the formation of EMSTO syst ems, such as the merging 
of tw o bound clusters (e.g., IMackev fc Brobv Nielsen! 
2007]) or the merger of a clust er and a giant molecular 
cloud (|Bekki fc Mackevl I2009D . fit with the observed 
properties of EMSTO systems. However, if the star 
formation has indeed progres sed unbroken over s evera l 
hu ndred Myr as sug g estedbv lRubele et all pOlOl I20T1 
or IGoudfrooii et al.l (|2011a|) . these scenarios are likely 
disfavored. 

If, as hypothesized, EMSTO clusters and Galactic 
globular clusters share a common formation process, a 
key prediction is that EMSTO clusters ought to har- 
bor similar star-to-star variations in the abundances of 
light elements as seen in the Galactic systems. A com- 
plicating factor is the comparatively high overall metal- 
licity of the EMSTO clusters ([M/H] -0.4), along 
with the possibility that their central gas resevoirs may 
have been augmented by an unknown amount of ac- 
creted pristine material (e.g., IConrov fc Spergel 120111) 
- so the scale of these variations is difficult to predict. 
Very few light element abundance measurements exist 
for EMSTO systems . The most extensive study is that 
bv iMucciarelli et al.1 ()2008h who targeted between 5 and 
11 stars in each of four intermediate- age LMC clusters, 
of which only one (NGC 1783) unambiguously features 
an EMSTO. For all four systems, Mucciarelli et al. as- 
sert that the star-to-star scatter in each of the ~ 20 ele- 
mental abundances they measure is negligible; in par- 
ticular, there are no O-Na or Mg-Al anti-correlations 
evident. However, the interpretation of these m e asure - 
men ts is somewhat c o ntrover sial - both iConrovl ([201 ll ) 
and IGoudfrooii et al.l (|2011bl) note that the star-to-star 
scatter in the listed abundances of sodium is larger than 
the observational uncertainties at a statistically signifi- 
cant level (up to ~ 4er) for some clusters. Even so, the 
sole EMSTO cluster in the study, NGC 1783, does not 
seem to exhibit elemental abundance patterns which are 
strikingly distinct from those of the other three systems. 

1.3. This work 

Additional high resolution spectral data, for a larger 
sample of stars in a larger ensemble of EMSTO clusters, 
are clearly required to assess the viability of the link be- 
tween these systems and the multi-population Galactic 
globular clusters. Such data possess an extra benefit be- 
yond exploring elemental abundance patterns. Because 
EMSTO clusters are diffuse, low-density stellar systems, 
their two-body relaxation times are typ ically comparable 
to, or longer than their ages (see e.g.. IGoudfrooii et al.l 
l2011b|) . Thus, unlike for Galactic globular clusters, the 
internal dynamics of EMSTO systems should still re- 
flect the conditions present early on in their lives - so 
that any signature imparted onto the cluster kinematics 
by the formation process should be both detectable and 
straightforward to interpret. 

This paper is the first in a series devoted to a de- 
tailed study of medium and high resolution spectra for 
21 giant stars and one planetary nebula in the most 



4 



Mackey et al. 



massive known EMSTO cluster, and one of the best 
studied photometrically - NGC 1846 in the LMC. Stel- 
lar populations in this cluster span the age range « 
1.6 — 1.9 Gyr, and its iron abundance is [Fe/H sa 
-0.4 ( e.g. jMackey fc Brobv Nielsenll2007t iMackev et ah! 
I2008aj) . The aim of our study of NGC 1846 is two-fold: 
first, to characterise the internal dynamics of the clus- 
ter and search for any signatures of its formation that 
might be present; and second, to place constraints on 
any star-to-star elemental abundance variations, espe- 
cially for light elements. Here, we present a detailed de- 
scription of our target selection and data analysis (Sec- 
tions [2] and 13]) , and focus on the cluster kinematics (Sec- 
tions [4] and [5]). Subsequent work will cover the elemental 
abundance analysis. 

2. OBSERVATIONS & DATA REDUCTION 
2.1. Data Acquisition 

We obtained spectra of stars in the vicinity of NGC 
1846 using the FLAMES instrument at the ESO Very 
Large Telescope (VLT) on Cerro Paranal, Chile, un- 
der programme 082 D-0387 (PI: Mackey). FLAMES 
(jPasquini et al.ll2002j ) is a fiber-fed multi-object spectro- 
graph mounted at the Naysmith A platform of the 8.2m 
Unit Telescope 2 (Kueyen). We employed the MEDUSA- 
GIRAFFE mode, allowing up to 132 stars to be targeted 
across the 25' diameter field of view in a single pointing. 

Our target selection is discussed in detail in Section 
12.21 We used just a single FLAMES fiber configuration 
but observed at three high-resolution GIRAFFE settings 
(HR11, HR13, and HR14B) and one low- resolution set- 
ting (LR02). The nominal wavelength coverage and spec- 
tral resolution for each of these set-ups is listed in Table 
[TJ This paper and the next in the series (Paper II) are 
dedicated to analysis of the observations obtained at the 
three HR settings; results from the blue LR set-up will 
be reported in a separate future work. 

We obtained our data in visitor mode on the three 
nights 2008 November 29 - December 01. Conditions 
were clear and stable, with seeing typically in the range 
0.5" — 1.0". On each night we observed a given HR set- 
ting at the beginning and end of the night, and reserved 
an hour either side of NGC 1846 crossing the meridian 
for observing the blue LR set-up with minimal differ- 
ential atmospheric refraction. For all four gratings we 
obtained 6 x 55 min exposures. Due to the faintness 
of our targets we switched the simultaneous calibration 
lamps off, instead bracketing each long science exposure 
with short calibration lamp exposures to ensure we could 
achieve an accurate wavelength solution. 

2.2. Target Selection and Photometry 

Targets for our FLAMES observations were drawn 
from the archival HST/ACS i maging of NGC 1846 de- 
scribed in the introduction (|Mackev fc Brobv Nielsen! 
120071: iGoudfrooii et al.ll200ll for the crowded central re- 
gions, and the Magellanic Clo uds Photometric Sur vey 
(MCPS) catalogue of the LMC (|Zaritskv et al.ll2004l ) for 
the surrounding field. 

For work such as this, one would ideally directly tar- 
get stars spread across the cluster EMSTO so that el- 
emental abundances and dynamical properties can be 
correlated directly against position on the CMD. Un- 



Table 1 

Nominal wavelength coverage and spectral resolution 
for our four GIRAFFE set-ups. 



Configuration 


Astart 




AA 


R 


Name 


(A) 


(A) 


(A) 




HR11 


5597 


5840 


243 


24 200 


HR13 


6120 


6405 


285 


22 500 


HR14B 


6383 


6626 


243 


28 800 


LR02 


3964 


4567 


603 


6 000 



fortunately, however, at V ~ 20.5, obtaining sufficient 
signal-to-noise for a high precision abundance analysis 
of NGC 1846 MSTO stars is beyond the capabilities of 
presently-available high-resolution spectrographs, even 
on the largest telescopes. Instead, we are forced to target 
the brightest giant stars in the cluste r. There is, however , 
no major disadvantage in doing so - IMilone et al.l (2009) 
demonstrated that the relative numbers of "younger" and 
"older" stars across the EMSTO are not too dissimilar (a 
ratio of rs 2:1). This means that, providing a sufficiently 
large ensemble of giant stars is observed, the full spread 
in population parameters (age and, if they exist, elemen- 
tal abundance variations) ought to be well sampled. 

At the time we were preparing our VLT observations, 
HST imaging of NGC 1846 was available from only two 
separate programs: 9891 (PI: Gilmore) and 10595 (PI: 
Goudfrooij), both of which utilised the Advanced Cam- 
era for Surveys (ACS) Wide Field Channel (WFC). The 
former is a "snapshot" program, where the imaging con- 
sists of just two exposures - 300s in the F555W filter 
and 200s in the F814W filter. The second program com- 
prises a more extensive set of observations - three expo- 
sures in each of the F435W, F555W, and F814W filters, 
where two of the exposures were long (340s each) and 
one short (90s, 40s, and 8s in the three filters, respec- 
tively). The two programs were observed with differ- 
ing orientations, meaning that their footprints only par- 
tially overlap on the sky. In the interests of covering as 
much of the area around NGC 1846 as possible with high 
quality ACS imaging and photometry, and hence max- 
imising the number of objects in our input catalogue for 
FLAMES target selection, we u sed the multidrizzle 
software ( Koe kemoer et all |2003|) to combine the com- 
plete set of F555W observations into a "master" refer- 
ence image. This image increases the sky coverage near 
NGC 1846 by ~ 30% over what would have been avail- 
able using only one of the HST programs^ 

We next used the DOLPHOT software (jDolp hin 2000), 
and in particular its ACS module, to photometer all the 
av ailable images. Details of this pro cedure may be found 
in IMackev fc Brobv Nielsen! (|2007l ) . Briefly, DOLPHOT 
performs point-spread function (PSF) fitting photometry 
using model PSFs especially tailored for the ACS cam- 
era. It works on images for which basic reduction steps 
have been applied (bias and dark current subtraction, 
and flat-field division) but which have not been distortion 
corrected (drizzled). The software can photometer mul- 
tiple images in multiple filters simultaneously, matching 
detections across images and deriving coordinates rela- 
tive to an input reference frame - in our case the master 
drizzled F555W image. Outp ut photometry i s on th e cal- 
ibrated VEGAMAG scale of iSirianni et all (f2005h and 
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Figure 1. HST/ACS color-magnitude diagrams for NGC 1846, in the (B — I, B) and [V — I, V) planes. In the upper two panels, the 
peculiar extended main-sequence turn-off morphology is clearly evident. The lower four panels show the regions on the red giant branch 
for our FLAMES target selection, bounded by dashed and solid lines, along with the objects we observed. Those targets marked with red 
dots are radial velocity members of the cluster (see Section [XT} ; those marked with blue crosses are not. 



has been corrected for charge-transfer efficiency degra- 
dation. Where possible, transformations are also made 
into the standard Johnson-Cousins system. To obtain a 
clean list of stellar detections with high quality photom- 
etry we filtered our dolphot measurements using the 
classification, sharpness and crowding parameters (see 
iMackev fc Brobv Nielsen! I2007t IMackev et al.ll2008al) . 

CMDs, in both the (B - I, B) and (V - I, V) planes 
are shown in Figure [TJ NGC 1846 is set against a rela- 
tively dense background of LMC field stars, so care was 
necessary in selecting our input catalogue for FLAMES. 
We used CMDs for stars within 30" of the cluster cen- 
ter to filter out most of the field contamination, allowing 
us to define regions encompassing the upper red-giant 
branch (RGB) and asymptotic giant branch (AGB) on 
both colour-magnitude planes. The colour and luminos- 
ity of a star on the RGB is sensitive to elemental abun- 
dances, especially that of iron, as well as to age. Given 
that the primary aim of our VLT program was to test 
for possible star-to-star abundance variations, and that 
we also strongly suspect the presence of an internal age 
spread in NGC 1846, we were particularly mindful not to 
define our target regions too restrictively. The red side 



of the NGC 1846 RGB is sharply defined on both CMD 
planes; however to the blue there is potential for sig- 
nificant overlap between the RGB and AGB sequences. 
Hence we set the blue side of our target region rather 
conservatively, with the result there are likely a few AGB 
stars in our FLAMES sample. As long as these are early- 
AGB objects (i.e., which have not yet undergone third 
dredge-up) the composition of the atmospheres of these 
stars should not have been altered substantially through 
the addition of newly-processed material. This issue is 
not relevant to the cluster kinematics, but is discussed in 
more detail in Paper II. 

We set an absolute red limit to our target region of V — 
I = 1.7 in order to exclude very cool giants for which it is 
more difficult to do a reliable abundance analysis. Note 
that this also excluded the most luminous (and evolved) 
AGB stars in the cluster. At the faint end we set a 
limit of V = 18.3, as spectra for stars fainter than this 
would not have sufficiently high signal-to-noise. On the 
(B — I, B) plane, the corresponding limits were B — I = 
3.25 and B = 19.35. 

Having defined our target regions we selected as our 
cluster input catalogue all stars across the full ACS field 
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Figure 2. Color-magnitude diagrams from the Magellanic Clouds Photometric Survey catalog (see Zaritsky ct al. 2004) for stars lying 
within the 25' diameter FLAMES field of view. Note that objects lying within the central w 2', where NGC 1846 is located, have been 
excised. The left panel shows the (J — K s , K s ) plane, where the photometry originally comes from the 2MASS catalog, while the right 
panel shows the (B — V, V) plane. Stars belonging to the LMC field for which we obtained FLAMES spectra are marked with red dots in 
both panels. Objects marked with blue crosses are our telluric correction stars (note that not all of these have infrared photometry). 



of view which lay within these regions on both colour- 
magnitude planes. As noted above, nearly a third of the 
field is not covered by F435W imaging - we filtered ob- 
jects in this area using only the (V — I, V) plane. The 
very outer parts of the ACS field were only sparsely pop- 
ulated with suitable stars, so we supplemented the in- 
put catalogue in this region with stars lying up to ~ 0.1 
mag to the red of our target region on the CMD. We 
also added to the catalogue a nearby planetary nebula 
likely belonging to NGC 1846, Mo- 17 (see Section ET2|) . 
We searched for, and removed from the catalogue any 
stars with neighbours within a radius of 2" that were 
sufficiently bright as to be likely to interfere with the 
spectrum of the target. 

We used ESO's Fibre Positioner Observation Support 
Software (FPOSS) to determine the optimal MEDUSA 
configuration for our FLAMES observations. A few trial 
runs using our ACS input catalogue demonstrated that 
we would be able to observe, at most, ~ 30 stars within 
« 2' of the center of NGC 1846 - the limiting factor being 
the minimum fiber-to-fiber separation of 11". Even al- 
lowing for ~ 15 sky fibers and ~ 5 fibers allocated to hot 
blue stars to allow correction of telluric absorption (see 
below), this left roughly 80 fibers unused. Rather than 
waste these we decided to allocate them to bright field 
RGB stars spread evenly over the non-cluster regions of 
the 25' diameter FLAMES field of view. 

We selected these objects from the Magellanic Clouds 
Photometric Survey (MCPS) catalogue of the LMC (see 
IZaritskv~ t al. 2004). This catalogue provides precise 
astrometry and UBVI photometry of stars in the cen- 
tral 64 deg 2 of the LMC, and, for some subset of suc- 
cessfully cross-matched stars, JHK S photometry from 
the Two Micron All Sky Survey (2MASS) catalogue 



(|Skrutskie et al.ll2006l ). There were a sufficiently large 
number of MCPS stars with 2MASS matches across our 
field of view that we were able to draw targets exclu- 
sively from this subset. The 2MASS measurements were 
particularly useful for distinguishing between field RGB 
and AGB stars, as shown in Figure [5] (left panel). We 
selected RGB stars brighter than V ~ 17 for our MCPS 
input catalogue (Figure [2] right panel). 

Available literature measurements suggested that the 
NGC 1846 rad i al ve lo city of V; vs g 240 kms -1 
HOlszewski et al.l 119911: IGrocholski et all 120061) would 
bring the forbidden [Oi] line at 6300.3 A into close prox- 
imity with telluric absorption features due to atmo- 
spheric 02- We planned to derive our oxygen abundances 
using primarily this line (see Paper II) ; in order to be able 
to account for the telluric absorption we required spec- 
tra of hot, preferably fast-rotating stars at high S/N. To 
this end, we selected a number of target stars from the 
MCPS catalogue that appeared to be bright members of 
the young blue LMC field main sequence (see Figure [5]) . 

To aid in the subtraction of atmospheric emission lines 
from our science spectra, we also decided to allocate ~ 15 
fibers to blank sky. We placed these such that for a given 
sky fiber there were no objects in the MCPS catalogue ly- 
ing within 5". Ideally these sky positions would be local 
to our highest priority targets (i.e., NGC 1846 members); 
however the very tightly packed nature of the fiber con- 
figuration in the middle of the FLAMES field meant that 
it was impossible to allocate any other fibers within ~ 5' 
of the cluster center. This ultimately led to some prob- 
lems achieving a high quality sky subtraction for many of 
our targets; however it was not difficult to accommodate 
these issues in our analysis ( Section 12.31) . 

Finally, we merged our ACS and MCPS catalogues 
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for input to FPOSS. The critical step was to transform 
the ACS coordinates onto the MCPS (FK5) astromet- 
ric frame. Although the relative astrometry of targets 
measured from the ACS imaging is extremely precise (of 
order milli-arcseconds), the HST FITS header informa- 
tion from which the absolute astrometry is derived can be 
in err or by up to several arcseconds (e.g. lAnderson et al.l 
120081) - more than enough for the MEDUSA fibres to 
miss these targets entirely. We cross-matched ~ 100 stars 
across the ACS master reference image with stars in the 
full MCPS LMC catalogue, excluding the crowded cen- 
tral cluster region, and used these to derive a suitable 
coordinate transformation onto the MCPS frame. 

Having successfully merged our two lists of potential 
targets, we used FPOSS to determine the optimal fiber 
configuration. In doing so, we assigned priorities to the 
various classes of target - for example, ACS stars received 
higher priority than MCPS stars and sky positions. The 
ACS stars themselves were graded in priority by lumi- 
nosity (the brighter the better) , distance from the center 
of the cluster (the smaller the better), and, weakly, by 
colour on the RGB (the redder the better). 

Our final configuration targeted 30 ACS stars (11 
within 50" of the cluster center), the planetary nebula 
Mo-17, 79 MCPS field RGB stars, 6 MCPS stars for cor- 
recting telluric absorption, and 16 blank sky positions. 
We observed our blue LR setting using FLAMES fiber 
positioner plate 1 and the three HR settings using plate 
2. On both plates, spectra from two fibers (both target- 
ing MCPS stars) fell off the edge of the CCDs and were 
lost. Additionally, two of the plate 2 fibers were broken 
leading to the loss of another MCPS star and one ACS 
target. Hence, the HR observations presented here re- 
sulted in spectra for 76 MCPS stars and 29 ACS stars, 
as well as the other targets and sky positions listed above. 

2.3. Data Reduction 

We used the ESO public GIRAFFE pipeline recipes 
v2.8.1, operating under the graphical front-end software 
GASGANO, to perform a basic reduction of all our sci- 
ence frames - that is, bias subtraction, fiber localization, 
optimal extraction of spectra, division by a normalized 
flat-field image (i.e., including correction for fibcr-to- fiber 
transmission differences), an initial wavelength calibra- 
tion, and rebinning to a uniform linearized dispersion 
scale (0.05 A per pixel for the spectra considered here). 
For each of the four instrumental set-ups, the science 
frames included the six on-target exposures as well as the 
short bracketing exposures with the calibration lamps on. 
For science frames belonging to a given set-up, the initial 
wavelength calibration was derived from an arc-lamp ex- 
posure taken at the beginning of the night on which the 
science frames were taken. After the basic data reduc- 
tion had been completed, we used the measured positions 
of emission lines in the bracketing calibration frames to 
check for small residual wavelength drifts in each of the 
on-target exposures individually. 

Next, we intended to perform a sky subtraction on each 
individual spectrum of an object, along with a correction 
for telluric absorption, before combining these spectra 
into a final product. We note that sky subtraction is 
only marginally important for the HR11 set-up, which 
covers a handful of weak emission lines; however the 
HR13 and HR14B set-ups both cover numerous bright 



emission lines. Similarly, telluric absorption is negligible 
(< 3%) for the HR11 and HR14B set-ups, but is notice- 
able for the HR13 set-up over the range 6275 - 6330 A 
where there are many lines of ~ 5 — 15% absorption due 
to atmospheric O2. 

Despite testing a variety of techniques we were unable 
to obtain a high-quality sky subtraction across the full 
wavelength coverage of either the HR13 or HR14B set- 
ups, especially for the likely members of NGC 1846. Us- 
ing the spectrum from the nearest sky fiber to a target 
did not work, nor did subtracting a combined spectrum 
derived from all 16 sky fibers in a given exposure. Us- 
ing the iraf task SKYTWEAK to vary the scaling of the 
sky spectra and apply small wavelength shifts improved 
the results marginally, but not to a satisfactory level. 
The origin of the problem is unclear, but it may well be 
linked to our inability to place sky fibers very locally to 
the NGC 1846 targets, as described in Section l2~2l 

To work around this issue we used the iraf task SCOM- 
BINE to merge all 96 sky spectra for a given set-up (16 
sky fibers from each of 6 exposures) into a high S/N 
"master" spectrum, and used this to generate a mask 
specifying all of the narrow wavelength intervals affected 
by sky lines. Across the wavelength range of the HR13 
set-up there were 33 such intervals in the mask, covering 
16.7% of any given spectrum, while for HR14B there were 
25 sky-line intervals excising 11.1% of the coverage and 
for HR11 there were 8 intervals masking 3.3%. When 
detecting and measuring stellar absorption features for 
our kinematic and chemical abundance analysis, we sim- 
ply ignored any lines lying within one resolution element 
(~ 0.3 A) of a masked sky region. As described in Pa- 
per II, this procedure did not reject any lines critical for 
examining the abundance of a given element, but did 
ensure that any lines impacted by sky emission did not 
adversely influence our measurements. 

We encountered a similar problem when attempting to 
correct for telluric absorption in the HR13 set-up. Four 
of our six bright blue targets turned out to be suitable 
stars for this type of correction; however despite trying 
a similar variety of techniques as for the sky subtrac- 
tion, we could not consistently achieve a high quality 
correction across the full affected wavelength range. Ulti- 
mately, we adpoted an identical approach to our sky-line 
solution - we normalised then combined the 24 telluric 
star exposures into a high S/N master spectrum and used 
this to generate a mask specifying the intervals affected 
by atmospheric absorption lines. There were 13 such in- 
tervals covering 10.1% of the HR13 wavelength range; 
we note some overlap between the telluric and sky-line 
masks. Again, during our kinematic and chemical abun- 
dance analysis we simply ignored any stellar lines lying 
within one resolution element of a masked telluric-line 
region. As described in more detail in Paper II, the only 
critical stellar line rejected by this procedure was the 
[Ol] line at 6300.3 A, which, 'for NGC 1846 targets, is 
redshifted onto a telluric feature near 6305.5 A. In order 
to use this line, we applied our master spectrum and the 
iraf task telluric to obtain a high-quality correction 
to only the 6305.5 A feature, and then used these spe- 
cially tailored spectra when deriving oxygen abundances. 

Finally, we corrected all on-target spectra to the helio- 
centric frame using the IRAF tasks RVCORRECT and DOP- 
COR. We then median-combined the six exposures of a 
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Figure 3. Representative spectral segments for three of our confirmed NGC 1846 members, as labelled. Each segment spans ss 15A. The 
left column contains a region from the HR11 coverage, the central column from the HR13 co vera ge, and the right column from the HR14B 
coverage. The spectra have been normalised following our DAOSPEC analysis (see Section 12.41 1 and corrected to the heliocentric frame; 
however no correction has been made for the individual radial velocities. Regions masked in our analysis due to the presence of sky lines 
are shaded yellow (note that in this example such regions appear in the HR14B spectra only). The top row of the plot corresponds to one 
of the brightest targets in the NGC 1846 sample, and the spectra are of commensurately high S/N. Similarly, the bottom row shows spectra 
for one of the faintest targets, and these are of correspondingly low S/N. The middle row contains data of intermediate S/N; two-thirds of 
our sample of confirmed NGC 1846 members have spectra of comparable or higher S/N. Note that even in the lowest S/N spectra, most 
absorption lines are still clearly visible. 



given target in a given set-up using the iraf task SCOM- 
bine. Best results (including cosmic ray removal) were 
achieved by scaling the input spectra to a common flux 
level and then, when combining them, weighting each by 
its median flux value and applying a sigma-clipping rejec- 
tion algorithm. Wavelength regions covered by the sky- 
and telluric-line masks were excluded when computing 
the scaling and weighting factors. 

Example spectra for three NGC 1846 stars spanning 
nearly the full range in brightness of our sample of con- 
firmed members (see Section 13- 1[) are shown in Fig. [3l 
We estimated the continuum signal-to- noise (S/N) for 
these objects with the iraf task SPLOT. The setting 
with the highest S/N per pixel is HR13; spectra for the 
HR11 and HR14B settings have S/N levels consistently 
~ 35 — 40% lower for the same exposure duration. Partly 
this is due to the slightly higher resolution of these two 
settings (the HR11 and HR14B settings span wavelength 
ranges s» 15% shorter than that for HR13 - see Table 
[1]) and partly it is due to their somewhat lower over- 
all efficiency The HR13 spectra span S/N ~ 30 — 70 
per pixel, while the HR11 and HR14B spectra have S/N 
in the range ~ 18 — 50 per pixel. For reference, the 
HR13 setting has « 5.6 pixels per resolution element at 
the mid-point of its coverage, whereas HR11 has « 4.7 
and HR14B 4.5. Overall, two thirds of our confirmed 
NGC 1846 stars have spectra with S/N comparable to, 
or greater than, the middle star plotted in Fig. [31 



2.4. Line Identification & Measurement 

We next applied the DAOSPEC software package to 
each of our 3 x (29 + 76) = 315 GIRAFFE HR spec- 
tra. DAOSPEC is an automated tool, optimised for 
spectra with R > 15 000, for identifying and measur- 
ing absorption lin e s. F ull details may be found in 
IStetson fe Pancinol (|2008f) . but briefly, it works by iter- 
atively finding lines in a given spectrum, fitting these 
with Gaussian profiles of fixed FWHM and subtracting 
them, and then using the residual spectrum to refine 
the continuum normalization and improve the line cen- 
troids and strengths. Once this process has converged, 
the measured lines are cross-correlated against a user- 
supplied line list to provide an estimate of the radial 
velocity of the target along with an identification for 
any lines that are successfully matched. The accuracy of 
DAOSPEC measurements, especially for equivalent widths, 
has been tested and verified on GIRAFF E HR data (e.g., 
iPompeia et~al"1l2008t ILetarte et al.ll2010D . 

Full details of our equivalent width measurements are 
provided in Paper II; here we are mainly interested in 
the radial velocities supplied by daospec. We assem- 
bled an input line list and atomic data for 18 neutral and 
12 singly ionised species (25 elements in total) over the 
wavelength int erval 5550— 6650 A using the VALD atomic 
line database (|Kupka et al.lll999f ). We included only the 
"strongest" lines for each species, where strength was 
approximately parametrized by the difference between 
the oscillator strength log(g/) and the excitation poten- 
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tial x (i- e -i by the sum log(g/) — x). For each species 
an empirical minimum limit for this value was defined 
by manually examining lines in our highest S/N spectra 
and identifying those with the smallest equivalent widths 
that could be reliably identified and measured (typically 
EW R3 15mA). Our final line list contained a total of 
1445 lines over the HOOA wavelength interval. Before 
input to DAOSPEC we applied our sky and telluric masks 
to cull affected lines from the list. 

When measuring the radial velocity of a target, 
DAOSPEC uses a sigma-clipping rejection algorithm to 
eliminate identified lines with discrepant velocities. In 
order to determine precise radial velocities for our stars 
we performed an initial run on all spectra with strin- 
gent 2cr rejection. This was important because at this 
stage we had not checked our line list for features that 
would be blended at the GIRAFFE HR resolutions. As 
described in Paper II, for the equivalent width measure- 
ments we subsequently performed a more relaxed 3.5cr 
rejection run and then determined the most appropriate 
lines for elemental abundance analysis by careful visual 
inspection of the spectra. However important additional 
guidance was provided by the identification of those lines 
initially rejected in the stringent 2er DAOSPEC run. 

For each of our targets we ultimately obtained three in- 
dependent radial velocity estimates from DAOSPEC, one 
for each HR setting. Typically, m 130 lines were identi- 
fied in each HR11 spectrum, w 100 lines in each HR13 
spectrum and 70 lines in each HR14B spectrum; however 
significant variation was seen as a function of S/N and 
the effective temperature of the target. In the majority of 
cases the line-by-line variance for each individual radial 
velocity estimate was in the range 0.5 < a < 1.5 kms , 
we conservatively adopted this quantity to represent the 
uncertainty on each such measurement. 

While checking our results we noticed the presence of 
small systematic offsets between velocities measured for 
a given star from different settings. The 2<7-clipped mean 
offsets were Vhri3 — Vhrii = 0.55 kms -1 (97 stars) and 
Vhri3 — Vhri4B = —0.36 kms -1 (94 stars). The origin 
of these offsets is not clear, but they were straightfor- 
ward to correct. For each setting we precisely measured 
the positions of all the atmospheric emission features in 
each of the 16 blank sky spectra. We then matched these 
m easurements against the wavelengths listed in the atlas 
of IQsterbrock et all (|1996[ ) and derived the mean (2.5<r- 
clipped) offset. We found that the HR11 spectra needed 
to be shifted by 0.39 kms' 1 , the HR13 spectra by —0.17 
kms -1 , and the HR14B spectra by —0.52 kms' 1 . No 
dependence on fiber number (i.e., the position of indi- 
vidual sky spectra on the FLAMES CCD) was evident. 
Applying these corrections almost completely removed 
the mean systematic offsets between velocities measured 
from the three different settings for a given target. We 
combined these corrected quantities in a weighted aver- 
age to obtain a final radial velocity for each star, along 
with its associated uncertainty. 

3. CLUSTER MEMBERSHIP 

3.1. Stars within the truncation radius 

The upper panel of Figure 0] shows radial velocity ver- 
sus projected distance from the center of NGC 1846 
for all 29 ACS stars and 76 MCPS stars. We adopted 
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Figure 4. Radial velocity versus distance from the center of NGC 
1846 for all of our FLAMES science targets. The upper panel spans 
the full radial coverage, and the full range of radial velocities. The 
vertical shaded area indicates the region lying inside the nominal 
boundary of the cluster, while the horizontal shaded region denotes 
our initial radial velocity cut of ±10 kms -1 about the mean ve- 
locity of targets lying within r% = 34.5". The lower panel zooms 
in on the intersection of these two regions. Here, the shaded areas 
indicate ±5 and ±3 times the velocity disper sion at given radius, 
calculated according to Eq. [4] (see also Figure li"0)l . In both panels 
the position of the planetary nebula Mo-17 is marked with a cross. 

a = 05:07:33.66, 5 = -67:27:40.7 for the cluster center, 
determined using our ACS photometry in the FK5 astro- 
metric frame. The tight grouping of objects to the lower 
left of the plot is indicative of the cluster. To determine 
membership we imposed a maximum allowed radius of 
161", correspondi ng to the truncat i on ra dius measured 
for NGC 1846 bv IGoudfrooii et al.l ([2009). We also cal- 
culated the mean velocity V of targets lying within the 
half-light radius of r^ = 34.5" obtained by the same au- 
thors, on the basis that these objects are the most likely 
to be cluster members, and then imposed a generous ra- 
dial velocity cut of ±10 kms" 1 about this valueQ- These 
criteria resulted in the exclusion of seven ACS stars, with 
one additional star lying on the boundary of the allowed 

3 We note that according to the 2010 update of the Harris ( 1996) 
online catalogue, all but the most massive Galactic globular clus- 
ters have central velocit y dispersions of only a few kms -1 (see also 
ILane et al.i r2009. 20103O. 
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Table 2 

Data on the radial velocity members of NGC 1846 studied in this work. 



Target 
Name 




Position (J2000.0) 
RA Dec 




Radius 

(") 


PA 

(°) 


B 


V 




Velocity 
(kms -1 ) 


Pmem 


ACS-001 


05 


07 


36 


.8 


-67 


27 


45. 


.9 


19.0 


105.8 


18 


.16 


16 


.56 


14 


.86 


236 


1 


± 


0. 


.5 


0. 


996 


ACS-013 


05 


07 


33 


.6 


-67 


26 


41 


.2 


59.5 


359.6 


18. 


.35 


16 


.99 


15 


.56 


238. 


.1 


± 


0. 


.5 


0. 


995 


ACS-017 


05 


07 


38. 


.4 


-67 


28 


11. 


.7 


41.2 


138.7 


18. 


.51 


17 


.18 


15. 


.73 


236. 


6 


± 


0. 


.4 


0. 


990 


ACS-025 


05 


07 


36. 


.2 


-67 


27 


58. 


.8 


23.2 


141.2 


18. 


.63 


17 


.28 


15. 


.93 


238. 


3 


± 





.5 


0. 


999 


ACS-030 


05 


07 


39 


.0 


-67 


28 


23 


.2 


52.5 


144.1 


18. 


.66 


17 


,13 


16 


.12 


240. 


7 


± 


0. 


.6 


0. 


985 


ACS-036 


05 


07 


30 


.4 


-67 


29 


35 


.7 


116.5 


189.4 


18. 


.86 


17 


.63 


16 


.38 


238. 


7 


± 


0. 


.5 


0. 


988 


ACS-043 


05 


07 


21 


.6 


-67 


27 


25. 


.3 


70.9 


282.5 


18. 


.84 


17 


.70 


16 


.45 


239. 


7 


± 


0. 


.5 


0. 


994 


ACS-046 


05 


07 


32 


.6 


—67 


27 


45. 


.5 


7.6 


230.5 


19. 


.11 


17 


.92 


16 


.68 


237. 


5 


± 


o 


.6 


0. 


997 


ACS-047 


05 


07 


34. 


.5 


-67 


28 


24. 


.2 


43.7 


173.7 


19. 


.02 


17 


.84 


16. 


.62 


239. 


2 


± 


0. 


.7 


0. 


997 


ACS-051 


05 


07 


36 


.6 


-67 


27 


33 


.5 


18.4 


66.8 


19. 


.10 


17 


.96 


16 


.71 


235. 


6 


± 


0. 


.8 


0. 


992 


ACS-053 


05 


07 


29 


.6 


-67 


26 


27. 


.4 


77.0 


342.3 


19. 


.09 


17 


.89 


16 


.65 


242. 


9 


± 


0. 


.8 


0. 


919 


ACS-059 


05 


07 


28. 


.6 


-67 


28 


11 


.8 


70.3 


204.3 


19. 


.15 


18. 


.01 


16. 


.85 


237. 


8 


± 


0. 


.8 


0. 


975 


ACS-066 


05 


07 


49. 


.6 


-67 


29 


01. 


.1 


121.9 


131.5 






18. 


.18 


17 


.00 


239. 


1 


± 


0. 


.8 


0. 


988 


ACS-080 


05 


07 


30. 


,3 


-67 


27 


11. 


.2 


35.3 


326.8 


18. 


.32 


17. 


.09 


15. 


.73 


243. 


9 


± 


0. 


.5 


0. 


974 


ACS-081 


05 


07 


30. 


.4 


-67 


28 


04. 


.2 


30.1 


218.8 


18. 


.45 


17 


.10 


15. 


.74 


239. 


2 


± 


0. 


.1 


0. 


998 


ACS-082 


05 


07 


30. 


.1 


-67 


27 


27. 


.4 


24.5 


303.1 


18. 


.47 


17. 


.11 


15. 


.81 


243. 


.4 


± 


0. 


.6 


0. 


984 


ACS-085 


05 


07 


30 


.0 


-67 


26 


42. 


.7 


61.7 


340.1 


18 


.53 


17. 


.15 


15. 


.81 


238. 


2 


± 


0. 


.6 


0. 


995 


ACS-090 


05 


07 


30. 


.2 


-67 


27 


46. 


.1 


20.8 


254.9 


18. 


.76 


17. 


.53 


16. 


.29 


238. 





± 





.6 


0. 


998 


ACS-092 


05 


07 


28. 


.7 


-67 


28 


24. 


.4 


52.2 


213.3 


18. 


.75 


17 
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Table 3 





Data on non- 


-members of NGC 1846. 




Target 


Position (J2000.0) 


Radius 


Velocity 


Name 


RA 


Dec 


(") 


(kms" 1 ) 


ACS-019 


05 07 14.5 


-67 28 16.4 


115.8 


278.7 ±0.5 


ACS-022 


05 07 20.0 


-67 29 36.9 


140.4 


309.9 ±0.7 


ACS-024 


05 07 21.1 


-67 27 01.3 


82.1 


249.1 ±0.5 


ACS-026 


05 07 23.6 


-67 26 50.6 


76.6 


262.8 ±0.5 


ACS-029 


05 07 14.7 


-67 29 15.8 


144.8 


279.5 ±0.5 


ACS-054 


05 07 26.3 


-67 26 56.8 


61.2 


292.0 ±0.8 


ACS-070 


05 07 10.5 


-67 28 31.1 


142.1 


281.0 ±0.8 


ACS-072 


05 07 14.4 


-67 27 48.1 


111.1 


291.0 ±0.9 



region of parameter space. The lower panel of Figure |4] 
shows a close-up of this region. Although the candidate 
star has a radial velocity only « 10 kms -1 higher than 
the systemic velocity of NGC 1846, this corresponds to 
more than five times the velocity dispersion at this radius 
within the cluster (see Section [4] below) and we therefore 
also excluded this object. 

Our 21 probable stellar members of NGC 1846 are 
listed in Table ED along with their ACS/WFC photome- 
try and measured radial velocities. Note, from the upper 
panel of Figure 21 that the velocity of the cluster overlaps 
significantly with the range measured for non-members 
(i.e., those objects outside the truncation radius). Thus, 
we cannot be certain that our assumed sample is entirely 
free of field interlopers. To assess the likelihood of this, 
we calculated an indicative membership probability for 
all 21 stars in our sample. For each object, we used the 
cluster surface density profile of lGoudfrooii et all (2009) 
(see the beginning of Section SI below) to estimate the 
likelihood that the star under consideration could be a 
non-member based on the relative densities of the cluster 
and field at the appropriate radius. We then counted how 
many of the_76 MCPS stars from outside r t lay within 
the interval V ± \ vi +0j|, where Vi and Oj are the velocity 
and associated uncertainty for the star under considera- 



tion, and used this information to estimate the likelihood 
that the star could be a non-member based on the devi- 
ation of its velocity from the cluster mean. We defined 
the final membership probability for the star, P mem , to 
be the complement of the product of this quantity and 
that estimated from the density profile. 

The values of P m em are listed in Table [2] Note that 
these are indicative lower limits only. They are subject to 
stochastic fluctuations due to the small number of MCPS 
objects lying comparably near to the cluster mean veloc- 
ity, and do not include any information about the pho- 
tometric selection criteria we employed in defining our 
initial sample of targets (Figure [T]) - which would effec- 
tively reduce the number of viable MCPS stars within 
the allowed velocity intervals by a factor > 2. 

The membership probabilities are very high for all ob- 
jects. The most likely non-member is ACS-053, which 
has a membership probability of ss 92%; for most others, 
Pmem ^ 99%. These simple calculations give us con- 
fidence that we have defined a clean set of NGC 1846 
members. We consider the remainder of our ACS and 
MCPS targets to belong to the LMC field (or, possibly, 
the Galactic foreground); further study of these objects 
is beyond the scope of the present work. 

In order to aid future observers of NGC 1846 we list, in 
Table [3l those stars lying within the truncation radius of 
161" which are not radial velocity members of the cluster. 
Members and non-members are also distinguished on the 
CMDs in Figure [1] 

3.2. The Planetary Nebula Mo- 17 

We briefly consider the planetary nebula Mo- 17 in 
more detail. This object was first catalogued bv iMorganl 
(1994), and, based on its proximit y to NGC 1846 (it lies 
at a radius of 83.5" 20.3 pc), iKontizas etHI (1996) 
suggested th at it might in fact bel ong to the cluster. The 
catalogue of iReid fc Parker! (|2006l ) lists a radial velocity 
of 241.8 ±20.0 kins" 1 which is consistent with literature 
estimates for the motion of NGC 1846. 
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Figure 5. HST ACS/WFC images of the NGC 1846 planetary 
nebula Mo-17 - from left to right, taken with the F435W, F555W, 
and F814W filters. Each thumbnail is 4" on a side. 



Figure |S] shows ACS/WFC thumbnails of Mo-17 in the 
F435W, F555W and F814W passbandfl The nebula is 
circular in appearance, with a diameter of approximately 
1.3" fa 0.3 pc. It is brightest in the F555W image, in 
which it also appears to have two distinct lobes, or per- 
haps an outer ring. The central star may be visible - 
there is a blend of two objects lying near the middle of 
the nebula. The fainter of these (the upper object in Fig. 
[5]) is much bluer than the brighter; however due to their 
close proximity, precise photometry is very difficult. 

As noted previously, we allocated a MEDUSA fiber to 
Mo-17. The fiber diameter of 1.2" corresponds well to the 
size of the nebula. Example spectra from the HR13 and 
HR14B settings are plotted in Fig. |6] Although there 
is sparse contamination from sky line residuals, in the 
HR13 setting [Oi] 6300 A emission is clearly visible, while 
in the HR14B spectrum there is strong emission from 
the [Nil] 6548A and 6583A lines and Ha. Also detected 
in the HR13 spectrum is a weak [Sm] 6312A line, and 
the [Oi] 6363 A line (not plotted). The HR11 setting (not 
shown) covers only one weak line, [Nil] 5755A, while the 
LR02 spectrum (also not shown) exhibits a variety of 
emission features including the Balmer lines FFy, US and 
He, several Hei lines, [Om] 4363A, and [Nelll] 3967A. 

In Fig. [6] is it quite evident that the strong [Oi] 6300A 
and [Nil] 6548A, 6583A lines are broad and double- 
peaked. The [Oi] 6363A line also follows this pattern. 
It is likely that this reflects the lobed structure seen in 
the F555W image in Fig. [5] The F555W filter spans 
the range ~ 4700 - 6100A, within which emission from 
[Om] 4959 A and 5007A will be dominant. The separa- 
tion of the peaks in the two [Oi] lines indicates outflow 
velocities of approximately ±20 kms , while the [Nil] 
peaks suggest ±10 kms -1 . 

The velocity structure present in these emission lines 
makes it difficult to obtain a precise estimate for the 
radial velocity of Mo-17. Notably however, the Ha line, 
while equivalently broad as the lines discussed above, has 
only a single peak and is highly symmetric. Ha falls be- 
tween the wavelength coverage of the ACS/WFC F555W 
and F814W filters. However, light in the F435W image 
of the nebula appears much more evenly distributed than 
in F555W; we know from our LR02 spectra that emis- 
sion m the range covered by this filter (~ 3700 - 4800A) 
comes predominantly from the Balmer series. Given this, 



4 Note that IShaw et al.l (2006) have previously reported on 
archival HST/ACS imaging of this object; however based on the 
coordinates and image they present, a nearby face-on background 
spiral galaxy was apparently misidentified as the planetary nebula. 
It is not clear whether their brief description of the appearance of 
Mo-17 corresponds to the correct object or this interloper. 
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Figure 6. Example HR13 and HR14B spectral segments for the 
NGC 1846 planetary nebula Mo-17. As in Fig. [3] these have been 
corrected to the heliocentric frame, but no correction for the radial 
velocity of the target has been made. Note that the intensity scale 
is in arbitrary units. We have subtracted the sky lines as accurately 
as possible across these small segments. As before, regions in the 
sky line mask are shaded yellow. Note that most lines have sub- 
tracted cleanly; however strong residuals are present for the [Oi] 
line near 6300 A. Regions masked due to the presence of telluric ab- 
sorption are shaded blue. The HR13 panel spans ~ 22A; nebular 
[Oi] emission is visible at 6305 A; note that this line is broadened 
and double-peaked. Weak [Sill] emission is also visible near 6317A. 
The HR14B panel spans a larger spectral range of ~ 45A. Double- 
peaked [Nil] emission lines are visible at either end of the segment, 
with a strong single-peaked Ha line near 6568A. This line was used 
to derive the velocity of Mo-17. 



and the circular symmetry of the nebula, we make the as- 
sumption that the peak of the Ha line reflects the overall 
radial motion of Mo-17, and derive a heliocentric velocity 
V r = 236.7 kms -1 . Estimating the mid-points of the two 
[Oi] and two [Nil] lines and taking the average of these 
measurements leads to a velocity within ss 1.2 kms -1 of 
this, and we thus adopt an overall uncertainty of ±1.5 
kms -1 on our result. As plotted in Figure [4j Mo-17 is 
very likely a member of NGC 1846; the indicative mem- 
bership probability, P mem ~ 97%. We include this object 
in Tabled 

4. CLUSTER KINEMATICS 

Our radial velocity measurements allowed us to inves- 
tigate the internal dynamics of NGC 1846. In what fol- 
lows we ma ke use of the clus t er str uctural parameters 
measured by Goudfrooii et al.l (|2009). These authors fit 
a [King (1962) model: 



n(r p ) = n 



1 



y/1 + (r P /r c ) 2 y/1 + (r t /r c f 



(1) 
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to a radial number-density profile constructed from their 
ACS/WFC imaging of the cluster. They found a core 
radius r c = 26.0" and a truncation radius r t = 161.2". 
They also measured r/, = 34.5" - assuming the cluster is 
not strongly segregated by luminosity (or mass), this is 
a good estimate of its projected half-light (or half-mass) 
radius. Adopting the canonical LMC distance modulus 
/i = 18.5, these values correspond to r c — 6.3 pc, = 8.4 
pc, and r t — 39.2 pc. 

4.1. Rotation 

We first checked for any signature of rotation within 
the cluster, following a com monly used proce dure (see 
e.g. JLane et^[2^0al2^10allbl ; lBellazzini et al.ll20lH and 
references therein). We calculated the position angle 
(PA, east of north) of each member in the plane of the sky 
with respect to the cluster center, then split the sample 
with the dividing line PA = 0° — 180° and calculated the 
difference in mean radial velocity between the two sub- 
groups of stars (AV r ). This process was then repeated 
with the position angle of the dividing line incremented 
by 20° each time. 

If coherent rotation is present in the cluster and does 
not lie entirely in the plane of the sky, this should mani- 
fest in the form of a sinusoidal pattern when AV r is plot- 
ted as a function of the position angle of the dividing line. 
The position angle at which the maximum amplitude of 
this sinusoid occurs corresponds to the projected axis of 
rotation, while the amplitude itself represents twice that 
of the mean rotation - i.e., A Iot = AV r /2. The observed 
(projected) amplitude A rot is a lower limit to the true 
amplitude of rotation A truc as there is a correction factor 
smi to consider, where i is the inclination of the rotation 
axis with respect to the plane of the sky (i = 90° repre- 
sents an "edge-on" cluster, where the observed amplitude 
of rotation would match the true amplitude, while i = 0° 
is a pole-on cluster where no rotation would be seen from 
radial velocity measurements irrespective of the true am- 
plitude). 

The upper panel of Figure [7] shows the results for our 
sample of 22 tracers in NGC 1846. We plot AV r /2 in 
order to more clearly elucidate the observed rotation am- 
plitude. To determine the point-to-point uncertainties, 
we used a bootstrapping method whereby we generated a 
large number of mock systems using our observed sample. 
Each mock system consisted of 22 stellar velocities cho- 
sen randomly from our 22 measurements, with repeated 
selection allowed. Every time we selected a velocity Vi 
we also generated a random Gaussian deviate of Vi (that 
is, a value from a Gaussian distribution with o~i - i.e., 
the measurement uncertainty in Vi) and added this to 
the velocity. We then repeated our measurement of AV r 
as a function of position angle for all mock systems, with 
the distribution of measurements at given PA indicative 
of the uncertainty at that point. This procedure allowed 
us to account naturally for both the (relatively small) 
size of our sample and the individual uncertainties on 
our measured radial velocities. 

Our results show a clear sinusoidal pattern, sugges- 
tive that systemic rotation is present in NGC 1846. The 
weighted best fit has A Iot = 1.1 kms" 1 and the axis of 
rotation lying along the line PA = 60° — 240°. Uncertain- 
ties derived directly from this fit likely under-estimate 
the true uncertainties, because the points are all corre- 
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Figure 7. Plots of AVr/2 as a function of the position angle of 
the dividing line for stars in NGC 1846. The upper panel shows 
the measurements and best-fitting model for the entire sample — 
note the clean sinusoidal shape of the rotation curve. The shaded 
region indicates the measured amplitude of rotation, .Arot — 

1.1 

kms . In the lower two panels we divide the sample in half using 
a radius r p = 55" and repeat the measurements. The inner half of 
the sample appears to rotate much more strongly than the outer 
half. For comparison to the upper panel, the shaded region again 
denotes the amplitude of rotation measured for the full sample. 



lated. iBellazzini et al.l ()2012f ) adopt ±0.5 kms" 1 and 
±30° as conservative la uncertainties on their measured 
rotation amplitudes and axis orientations, based on ex- 
perimentation with sub-samples of their observed stars 
and comparisons with external samples. We are not in a 
position to conduct similar tests; instead we investigated 
this problem using another large set of mock systems. 

The question is how strongly we can constrain the am- 
plitude and orientation of rotation in the cluster given the 
sample of stars we have measured. To this end each mock 
system consisted of 22 stars at the same radius and PA as 
in the real sample. We assumed simple cylindrical rota- 
tion with amplitude and orientation as measured above 
- this defined a velocity for each star based on its PA, to 
which we added a random measurement uncertainty as 
previously, and a second random component defined by 
the velocity dispersion profile of NGC 1846 (see Section 
14.21 below). We then followed the same procedure out- 
lined above to measure the rotation amplitude and orien- 
tation, and examined the distributions of these recovered 
quantities over all mock systems to determine our uncer- 
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tainties. We found the measured amplitude of rotation 
to be uncertain at the ±0.4 kms -1 level, and the orien- 
tation of the rotation axis to be uncert ain at the ±20° 
level - similar to the values adopted by iBcllazzini et al.1 
(l20H . 

We next decided to test whether there was any change 
in the degree of rotation with radius within the cluster. 
To do this we split our sample in half using a dividing 
radius r p = 55", which corresponds to roughly 1.5 times 
the half-light radius, and repeated our measurements us- 
ing only stars within, and then outwith, this radius. The 
results are shown in the lower panel of Figure [7] Al- 
though the samples are by now rather small, there is an 
indication that the rotation amplitude for the inner sam- 
ple is significantly larger than that for the outer sample 
- the best fits have amplitudes of 1.8 ± 0.6 kms -1 and 
0.8±0.4 kms -1 respectively. The PA of the rotation axis 
changes by only a small amount for the inner sample, to 
47 ± 25°, but increases substantially for the outer sam- 
ple, to 109 ± 35°. The quoted uncertainties were again 
d erived using larg e sets of mock systems. 

iBellazzini et al.l ((2012) note that variation of rotation 
amplitude with radius is not unusual in Galactic glob- 
ular clusters, such that the rotation amplitude deter- 
mined using stars at all radii typically under-estimates 
the maximum rotation amplitude in the cluster. Their 
experiments indicate that AV r = 2A TOt , where A rot is 
calculated (as we have done here) using stars at all radii, 
is likely to be a better estimate of the maximum rota- 
tion amplitude. This is consistent with the lower panel 
of Figure [71 which suggests that the maximum projected 
rotation amplitude for NGC 1846 could be as high as « 2 
kms -1 . 

As a qualitative check on our detection of rotation in 
NGC 1846 we used the best-fit solution derived from the 
full sample to construct a rotation curve. The top panel 
of Figure [5] shows a tangent plane projection of the po- 
sitions and velocities of our measured members within 
the cluster. The preferred axis of rotation, lying along 
the line PA = 60° - 240° ± 20°, is indicated. In the 
lower panel we plot, for each star, the offset from the 
cluster systemic velocity (derived in Section 14.21 below) 
as a function of the perpendicular distance from the rota- 
tion axis (Apa ). Bearing in mind that almost all of the 
point-to-point scatter can be accounted for by the cluster 
velocity dispersion (which is of order 2 kms -1 - again, 
see Section 14.21 below), these plots clearly support the 
detection of systemic rotation in NGC 1846. As an in- 
dica tive measure, we fit a rotation curve of the form (as in 
g..lLvnden-Bel]ill967l:lGottlll973l:rHenault-Brunet et al.l 



I20'l2alk 



2A rot 

Vrot = X 



x 



PA 



fpcak 1 + (ApA /V peak ) 2 



(2) 



where r pca k is the projected radius at which the peak 
amplitude of rotation occurs. The best fitting model 
has A rot = 1.8 ± 0.7 kms -1 and r pca k = 4.4 ± 3.1 pc. 
Although not particularly strongly constrained, these re- 
sults are entirely consistent with those from our initial 
analysis above. 

4.2. Mean velocity and velocity dispersion 

Next, we corrected the velocity of each of our tracers 
for the systemic rotation, and calculated the global mean 
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Figure 8. Additional evidence for systemic rotation in NGC 1846. 
The upper panel shows a tangent plane projection of the positions 
and velocities of the measured stars within the cluster. Red crossed 
points (alternatively blue dotted points) indicate stars with veloci- 
ties greater (less) than the global mean velocity. The point sizes are 
proportional to the deviation of the velocity from the mean. The 
preferred axis of rotation, oriented at a position angle 60 ± 20° 
east of north, is marked. In the lower panel we plot a rotation 
curve for NGC 1846, where -XpA is the perpendicular distance 
of from the rotation axis and we have assumed an LMC distance 
modulus fj, = 18.5. The best fitting model, of the form of Equa- 
tion[2]is marked, along with the curves represented by ±1<t on the 
maximum amplitude (see text). 



velocity and velocity dispersion, along with the velocity 
dispersion as a function of distance from the cluster cen- 
ter. To make the rotation corrections, we used the best 
fit to the full-sample rotation curve (i.e., the top panel of 
Figure [7J . Although, as discussed, we suspect that the 
degree of rotation does vary as a function of radius within 
the cluster, our sample of kinematic probes is not sufH- 
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Figure 9. Relative likelihood distributions for the systemic veloc- 
ity Vsys and the mean velocity dispersion a c i, where in each case we 
have marginalized with respect to the other variable. The dashed 
horizontal lines indicate the levels corresponding to la, 2a, and 3<r 
uncertainties. 

ciently large to accurately measure this variation; the 
use of the mean rotation curve represents an adequate 
compromise. In any case, we note that even the extreme 
case of applying no correction for rotation changes the 
global mean velocity by less than the uncertainty on the 
measurement, and the velocity dispersion results by less 
than ~ 15%. 

To calculate the mean radial velocity and velocity dis- 
persion for NGC 1846 we used a maximum likelihoo d 
technique following that defined by I Walker et al.l {2006). 
We assume that the measured velocities for our stars, de- 
noted {vx, . . . , vn}, are normally distributed about the 
systemic velocity V sys according to the associated mea- 
surement uncertainties Oi and the intrinsic cluster veloc- 
ity dispersion a c i . We can obtain numerical estimates for 
the quantities V sys and a c i by maximising the logarithm 
of the joint probability function for {v\, . . . , vn} - i-e., 

1 N 1 N ( — V \ 2 1 

Hp) = -\ E^?+^)-~ E % -\*H**). 

Z i=l Z i=l ^ +a cl> 1 

(3) 

Figure [9] shows the relative likelihood of both V sys and 
G c i , where in each case we have marginalized with respect 
to the other variable. To determine la uncertainties we 
calculated the parameter limits for the region containing 
the central 68.3% of the distribution function. 

We find a systemic velocity for NGC 1846 of V sys = 
239.1 ± 0.4 kms" 1 . This is con sistent with the value o f 
240 ± 10 kms" 1 determined by IQlszewski et al.l (|1991[) . 
but slightly large r than that of 235.2 ± 0.9 kms" 1 from 
IGrocholski eTafl J2006). For the mean velocity disper- 
sion we find a c i = 1.81^g'29 kms" 1 . A value of zero is 
excluded at ^> 5a, indicating that we have resolved the 
intrinsic cluster dispersion. 

The kinematic properties of globular clusters are com- 
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Figure 10. Our indicative three-point velocity dispersion pro- 
file for NGC 1846, together with the best-fitting Plummer model 
described in the text. The shaded region corresponds to the la 
uncertainties on the central dispersion ao- 

monly framed in terms of the central velocity dispersion, 
<7q. Given that we have members extending to ~ 4r/j it is 
likely that our measurement of a c i is an under-estimate 
of this quantity. To assess this we recalculated the veloc- 
ity dispersion as a function of cluster radius. Because of 
our relatively small sample size we are restricted to just 
three radial bins, corresponding roughly to stars inside 
r/j, stars between and twice r^, and stars outside 2r^. 
For each bin we utilized Eq. [3] as before, but with fixed 
y sys = 239.1 kms" 1 . 

Our results are shown in Figure [T0l A marginally sig- 
nificant decrease in the velocity dispersion with radius 
is evident. Even in the outermost bin the dispersion 
is resolved (i.e., a value of zero is excluded at a « 3u 
level). As an indicative meas ure we assume N GC 1846 
is isotropic and fit a pro iected iPlummen (|1911[ ) model to 
our three-point dispersion profile: 

a\r v ) = ^^ = . (4) 

Here a is a sca l e radi us which, for the family of pro- 
jected |Plummei] (|1911D models, is equal to the half-light 
(half-mass) radius Th if the mass-to-light ratio M/L is 
constant within the cluster. Note that we are not assert- 
ing here that a Plummer model is the most appropriate 
model to describe the internal kinematics of NGC 1846 
- with our poor radial resolution we are not in a posi- 
tion to undertake such an analysis. Rather, we use the 
Plummer model as a convenient parametrization to in- 
vestigate how (To relates to both our measurement of a c i 
and our measure of a within ss r^, in a system with a 
constant-density corc0. 

In principle both ao and a are free parameters in Eq. |4j 
The degree by which the scale radius a differs from ?'/, as 

5 Most globular clusters are found to conform closely to models 
in which the density of stars within the cluster core is approxi- 
mately constant. The family of [King (1962, 1966]) mo dels are the 
most pro minent exampl es, but others such as those of El son et alj 
(1987) or Wilson (1975) are seen to provide s uperior fits in some 
c ases (e.g., [M cLaughlin & van dcr Marcl 2009). The radial profile 
of IGoudfrooii et all l|2009l l for NGC 1846 clearly shows this cluster 
to be well described by models possessing constant-density cores. 
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determined from a surface-density or surface-brightness 
profile indicates how M/L may va ry with radius in the 
cluster (see e.g.. lLane et al.ll2010b[ ). Here, however, we 
again do not have a sufficiently large sample of kinematic 
tracers, especially at inner radii, to put adequate con- 
straints on a. We therefore make the assumption that 
mass follows light (i.e., a = — 34.5") and fit for cjq. 
The results of this process are seen in Figure [K3 We 
find (To = 2.521°'^ kms -1 , only slightly larger than the 
value measured for our centralmost bin, but considerably 
in excess of our global dispersion. 

4.3. Cluster mass & luminosity 

Having estimated the central velocity dispersion of 
NGC 1846 we can in turn derive an estimate o f its mass, 
since , for an isotropic Plummer model (see e.g.. lDeionghel 
[1981 : 



M= 64oo| 
3ttG 



(5) 



Substituting in our assumed value for a — 8.4 pc and our 
estimate for cto, we obtain M = (8. 4l^) x 10 4 M Q . 

We would also like to obtain an estimate of the global 
mass-to-light ratio for NGC 1846. To do this we first 
need an expression for the cluster luminosity, which we 
derive using Eq. Q] integrated with respect to 2irr p dr p : 



L{r p ) 



J V,0 



In (a) 



y/a— 1 

VP 



1 



(6) 



where 



a = 1 + {r p /r c f and = 1 + (r t /r c ) 2 . 

Here, £y,o is the V-band surface brightness correspond- 
ing to no in Eq. [T] Not e that in directly scalin g the 
number-density profile of iGoudfrooii et all (|2009) to a 
surface-brightness profile we assume that NGC 1846 is 
not strongly mass segregated. That it is a very dif- 
fuse cluster for i ts age (see e.g., iMackev et al.l [2008b; 
iKeller et al.ll20lil ) is consistent with this assertion. Note 
that we also made this assumption previously when set- 
ting a — rh in Eq. [4] in order to estimate do- 

To obtain a measurement of the total cluster y-band 
luminosity Ly we set r p — r t in Eq. [5] and estimated 
the central surface brightness Sy(0) directly from our 
master ACS/WFC F555W reference image. To mitigate 
the effects of any mis-centering, and random fluctua- 
tions due to the brightest stars, we measured the flux 
within a relatively large aperture of radius r p = 10". 
This is acceptable because the core radius of NGC 1846 
is even larger still. We obtained IV(0) = 19.6 ± 0.1 
mag arcsec" 2 , which, assuming the LMC distance mod- 
ulus is 18.5 and the foreground reddening E(B — V) = 
0.08, corresponds to Sy(0) = 670 + 62 L© pc~ 2 . How- 
ever, we must also take into account that NGC 1846 
is set against a moderately dense LMC field - we mea- 
sure the local surface brightness to be XVbfcg = 22.8 
mag arcsec~ 2 = 35 Lq pc -2 by using regions on the 
ACS/WFC image beyond r t . Thus, for the cluster only, 
£ y (0) = 635 ± 62 L Q pc~ 2 and, substituting into Eq. [IJ 
we find Ey-,0 = L41Sy(0) = 897+ 87 L Q pc" 2 . Then, 
from Eq. L v = (1.44 ± 0.14) x 10 5 L . 

This leads to a global mass-to-light ratio M/Ly = 

°- 59 iC:io for NGC 1846 - To p lace this value in a use- 



ful context, we compare it to the predictions of two 
well-known pop ulation synthesis codes - the Pe gase 
(v2.0) models of iFioc fc Rocca-Volmerange] (2997), and 
the Galaxev models of iBruzual &: Chariot! (|2003f ) - for 
a cluster of the age and metal abundance of NGC 1846. 
For both sets of models we assume Z — 0.008 and a single 
instantaneous burst of star formation. We further adopt 
the default initial stellar mass ranges - 0.1 — 120 M Q for 
Pegase and 0.1 - 100 M Q for Galaxev - and the clos- 
est av ailable initial mass function (IMF) to the iKroupal 
(|2001| ) "universal" IMF. For the Pegas e mode ls this cor- 
responds to the IMF of lKroupaet al.l (11993f) while fo r 
the Galaxev models it is the IMF of lChabrierl (|2003|l 
The Pegase models allow us to make the assumption 
that all white dwarfs are retained in the cluster, while 
all neutron stars and black holes are expelled upon for- 
mation. For the Galaxev models we necessarily assume 
all stellar remnants remain in the cluster; the predicted 
mass-to-light ratios reported below would drop by a few 
hundredths if the contributions from neutron stars and 
black holes were excluded. 

The results are as follows. The Pegase models predict 
M/L v = 0.69 at an age of 1.6 Gyr and M/L v = 0.88 
at 2.0 Gyr, while the Galaxev models predict M/L v — 
0.53 and 0.70 at these two ages. The differences between 
the two sets of results seem to be mainly due to the 
assumed IMFs and initial stellar mass ranges. Taken 
at face value and within the precision of our present 
measurements, the global mass-to-light ratio we have 
observed for NGC 1846 agrees acceptably with expec- 
tations derived purely from analysis of its constituent 
stellar populations. Note, however, that the systemic 
rotation we have detected could well play an important 
enough role in the internal kinematics of the cluster that 
the use of Eq. [5] when estimating M/Ly may not be 
appropriate (see Sections 14.51 and 15]) . 

4.4. Relaxation time 

Finally, we take the parameters measured above and 
use them to derive the central a nd half-mas s relax ation 
times in NGC 1846. Following iDiorgovskil (| 19930 . the 
central relaxation time is given by: 



t rc = 8.338 x 10 b yr 



1/2 ■j 
PO r c 

mln A 



(?) 



where po is the central mass density in the cluster, m is 
the mean stellar mass, and the units of mass and distance 
are solar masses and parsecs, respectively. The quantity 
In A is the Coulomb logarithm where A 0.4 A if N is 
the total number of stars in the system. 

Simpl e population synth esis models (such as SIM- 
CLUST, iDeveikis et al.l I2008D show that in a cluster of 
age w 1.75 Gyr with a IKroupal (|2001l ) mass function, 
the mean stellar mass in sa 0.5 M . Here we again as- 
sume that NGC 1846 is not strongly mass segregated so 
that m does not vary significantly throughout the system. 
The number of stars in NGC 1846 is N = M/fh 2M. 
To derive p n we used the two procedures outlined by 
iDjorgovskil |l993l ) to convert our observed, extinction- 
corrected central surface-brightness IV (0) to a central 
luminosity density jo = 71 + 5 L Q pc~ 3 , and then 
multiplied by our derived mass-to-light ratio so that 
po = 42lg° M pc~ 3 . Substituting all these values into 



16 



Mackey et al. 



E q . m 

we obtain a central relaxation time t rc — 2.4^Q 2 

Gyr. 

Th e half-mass, or median, re laxation time is given by 
(e.g., iBinnev fe Tremaindll987f ): 

t rh = 2.055 x 10 6 yr x _ - h k . (8) 

mm A 

As before, the units of mass and distance in this equation 
are solar masses and parsecs, respectively. Application of 
the appropriate values for NGC 1846 yields t r h = 2.6+ ; 2 
Gyr. This is only marginally longer than the central re- 
laxation time, consistent with the observation that NGC 
1846 is a diffuse cluster (i.e., the core radius r c is com- 
parable to r/j). 

4.5. Robustness of the kinematic analysis 

The magnitude of the rotation we infer in NGC 1846 
is roughly comparable to the mean velocity dispersion in 
the system, which is quite unusual for a stellar cluster 
(see Section [5j below) . Because of the potential signifi- 
cance of this observation, it is important to thoroughly 
assess the robustness of our kinematic analysis. There 
are two main reasons why it is necessary to undertake 
such a check. First, our sample of dynamical probes is 
comparatively small and so it is not inconceivable that, 
for example, a chance arrangement of velocity and posi- 
tion angle among a few stars could lead to a false rota- 
tion detection. Second, NGC 1846 has a non-neglibigle 
binary star fraction (note the strong binary star main 
sequences seen in the two CMDs in Figure [T]) and so for 
a subset of stars in our ensemble the observed velocities 
likely possess extra components due to stellar compan- 
ions rather than constituting a pure representation of the 
cluster kinematics. 

To investigate these issues we developed 17 sets of 10 5 
random realizations of our measured sample, and calcu- 
lated for each set (i) how frequently we detected a rota- 
tion signal comparable to that observed for NGC 1846; 
and (ii) the mean ensemble velocity dispersion. We char- 
acterized a given set of models according to the assumed 
properties of its binary star population. This analy- 
sis is quite similar, a t least in its principle, to that of 
IMcConnachie fe Cotel (l2Moh . 

Each realization consisted of 22 targets at identical 
radii to our real sample, but at random position angles 
and with random velocities generated according to a well- 
defined set of rules. It could be argued that the radius 
assigned to each star should also be randomly selected, 
according to the density profile of NGC 1846. However, 
our measured sample is not solely defined by this den- 
sity profile - rather it is the profile convolved with some 
complex and essentially unknown selection function for 
input to FLAMES (see Section [272]) . In the absence of 
this information, we felt the most sensible option was to 
maintain the radial distribution of the observed sample 
for each mock system. 

For each target we first assigned a base velocity by se- 
lecting randomly from a Gaussian distribution of width 
specified by Equation@]with a = = 34.5" and <Jo = 2.5 
kms -1 , evaluated at the appropriate radius. To this we 
added a random Gaussian deviate to repre sent the obser- 
vational uncertainty, exactly as in Section 14.11 Next, we 
randomly assigned the star to be an unresolved binary, 



or not, according to the assumed binary fraction for the 
overall set to which the mock realization belonged. If a 
star was selected to be a binary we randomly generated 
a set of parameters to characterize the system, again ac- 
cording to the properties of the overall set, leading to 
an additional component to add to the base velocity. If, 
however, a star was not selected to be a binary, no fur- 
ther modification of its base velocity was made. Note 
that the planetary nebula Mo- 17 was never selected to 
be a binary (since our velocity for this object comes from 
the nebula itself rather than the central star), but was 
otherwise treated identically to the other 21 stars in the 
sample. 

Binary systems are characterized by the masses of the 
two components, mi and mi (where the mass of the pri- 
mary mi > m,2), the orbital period P, and the orbital 
eccentricity e. It is convenient to define the mass ratio 
q = ma/mi. The semi-major axis ai of the orbit of the 
secondary about the barycenter may be calculated using 
Kepler's third law, and the semi-major axis of the pri- 
mary a\ = qa-i. To place a given binary system into the 
observational plane requires several additional parame- 
ters - the inclination i of the orbit to the line of sight, 
the argument (or longitude) or periastron, oj, and the 
orbital phase 9 at which the observation was made. The 
radial velocity of the primary then varies as: 

2irai sini , , 
Vb,i = — 1= [cos(0 + u) + e cos(w)] . (9) 
Py 1 — e 2 

We defined each set of random realizations according 
to the distributions from which the mass ratio, period, 
and eccentricity of each binary were drawn. Mass ra- 
tios were selected from eit her a uniform distribution o r 
the normal distribution of iDuquennov fe Mavorl (|1991l ). 
while periods were chosen according to either a uniform 
distribution in logP or th e log-normal distribution of 
IDuquennov fe May or (1991). Orbital eccentricities were 
selected from one of four distributions: circular orbits 
only (e = 0), a uniform distribution, the normal dis- 
tribution of IDuquennov fe Mavorl ([19911 ). or the thermal 
distribution of Heggie (1975). The different combinations 
of these distributions defined 16 sets of models. We also 
computed a seventeenth set with no binaries for compar- 
ison purposes. 

Since all our observed targets are upper RGB stars 
(with a few possible AGB stars), for every mock bi- 
nary system we set the primary mass mi = 1.625Af Q , 
based on the typical masses of upp er RGB stars in the 
best-fitting Da rtmouth isochrones dpotter et ail 120081) 
for NGC 1846 (|Goudfrooii et al.ll2009fl . The binary frac- 
tion for a gi ven set of models w as defined using the ob- 
servations of lMilone et a l. (2009) extrapolated according 
to the assumed distribution of mass ratios. Milone et al. 
estimated the binary fractions for 12 intermediate- age 
LMC clusters from HST/ACS imaging by counting stars 
above the main sequence. Due to blends and photometric 
uncertainties they were only able to clearly identify bina- 
ries with q > 0.6 or 0.7, depending on the quality of the 
imaging. They found typical binary fractions fb ~ 0.14 
for q > 0.6, or /& ~ 0.09 for q > 0.7. For a a uniform 
distribution in q these corre spond to a total binary frac - 
tion ft, ss 0.33, while for the IDuquennov fe Mavorl (|199lD 
normal distribution they correspond to fb ~ 0.55. We 
adopted these two values for our various sets of models, 
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Table 4 

Results from the mock realizations of our measured ensemble, including the effects of unresolved binary stars. 



Set 


Assumed distributions 


Peak > 1 . 1 km s 


- 1 (96) 


Pcak> 0.7 km s" 1 (%) 


Mean kinematics (kms 1 ) 


Number 


q logP e f b 


A ro t + rms 


+R 2 


A rot + rms +R 2 


V ByB °V ayB &cl &o cl 



01 








0.00 


2.01 


0.23 


0.39 


8.88 


1.31 


0.97 


239.10 


0.41 


1.79 


0.32 


02 


Unif. 


Unif. 


Circ. 


0.33 


6.28 


0.33 


0.72 


18.15 


1.35 


1.30 


239.11 


0.51 


2.18 


0.42 


03 


Unif. 


Unif. 


Unif. 


0.33 


5.98 


0.34 


0.74 


17.80 


1.34 


1.23 


239.09 


0.50 


2.18 


0.42 


04 


Unif. 


Unif. 


Ther. 


0.33 


5.93 


0.32 


0.75 


17.66 


1.32 


1.26 


239.11 


0.50 


2.16 


0.43 


05 


Unif. 


Unif. 


Norm. 


0.33 


6.20 


0.33 


0.73 


17.96 


1.26 


1.23 


239.09 


0.50 


2.18 


0.42 


0(5 


Unif. 


Norm. 


Circ. 


0.33 


6.01 


0.33 


0.69 


18.10 


1.35 


1.26 


239.11 


0.50 


2.16 


0.42 


07 


Unif. 


Norm. 


Unif. 


0.33 


5.98 


0.31 


0.67 


17.58 


1.38 


1.26 


239.10 


0.50 


2.17 


0.42 


08 


Unif. 


Norm. 


Ther. 


0.33 


5.86 


0.38 


0.76 


17.44 


1.30 


1.25 


239.10 


0.50 


2.16 


0.41 


09 


Unif. 


Norm. 


Norm. 


0.33 


5.76 


0.31 


0.75 


17.61 


1.29 


1.25 


239.10 


0.50 


2.17 


0.42 


10 


Norm. 


Unif. 


Circ. 


0.55 


8.21 


0.32 


0.90 


21.73 


1.21 


1.40 


239.10 


0.53 


2.33 


0.43 


11 


Norm. 


Unif. 


Unif. 


0.55 


8.24 


0.33 


0.86 


21.68 


1.25 


1.34 


239.09 


0.53 


2.33 


0.43 


12 


Norm. 


Unif. 


Ther. 


0.55 


8.07 


0.33 


0.89 


21.58 


1.21 


1.38 


239.10 


0.54 


2.33 


0.43 


13 


Norm. 


Unif. 


Norm. 


0.55 


8.23 


0.38 


0.88 


21.72 


1.25 


1.33 


239.10 


0.54 


2.32 


0.42 


14 


Norm. 


Norm. 


Circ. 


0.55 


7.96 


0.34 


0.81 


21.14 


1.24 


1.30 


239.09 


0.54 


2.30 


0.42 


15 


Norm. 


Norm. 


Unif. 


0.55 


7.78 


0.32 


0.89 


21.19 


1.26 


1.40 


239.10 


0.52 


2.29 


0.43 


16 


Norm. 


Norm. 


Ther. 


0.55 


7.80 


0.36 


0.90 


20.96 


1.28 


1.39 


239.09 


0.53 


2.31 


0.42 


17 


Norm. 


Norm. 


Norm. 


0.55 


7.79 


0.29 


0.80 


20.93 


1.22 


1.31 


239.10 


0.53 


2.29 


0.42 



depending on the assumed distribution in q. 

Having selected a given star to be a binary, our pro- 
cedure for generating a radial velocity was as follows. 
We first randomly selected a mass ratio, period, and ec- 
centricity from the relevant distributions. We defined a 
lower bound to allowed values of q by noting that the 
mass limit for hydrogen-burning, mi « O.O8M0, corre- 
sponds to q w 0.05 in our case. We further defined limits 
to the a llowed values of P following the procedure out- 
lined bv iMcConnachie fc CoTel ([2010). Briefly, the lower 
limit is set by the minimum orbital separation before 
the onset of mass-transfer, while the upper limit is set 
by the orbital separation corresponding to the bound- 
ary between "hard" and "soft" binaries in NGC 1846 - 
the latter of which are efficiently destroyed by three- or 
four-body interactions within the cluster. Next, we de- 
termined the orientation of the binary by selecting the 
angles i and lo from uniform distributions. Finally, we 
selected a number in the range 0—1, also according to 
a uniform distribution, to represent the elapsed fraction 
of the binary's orbital period since periastron, and con- 
verted this into the phase 9 by numerically solving Ke- 
pler's equation. With all the necessary parameters in 
hand, we determined the radial velocity according to Eq. 
[9l To mimic our original identification of cluster mem- 
bers (Section 13. ip . if the resulting velocity of the star 
was separated from the systemic cluster velocity by more 
than five times the velocity dispersion at that radius, we 
returned and regenerated all its binary parameters. 

The results for all 17 sets of realizations are presented 
in Table |4j To illustrate our analysis procedure we con- 
sider set 1 - the control set with zero binary fraction. We 
quantified whether the inferred rotation of a mock sys- 
tem matched that observed for the real cluster, by using 
the amplitude of rotation and two quality of fit param- 
eters. As described in Section 14.11 the rotation velocity 
implied by the NGC 1846 curve shown in the top panel 
of Figure [7] is 1.1 ± 0.4 kms^ 1 . The curve is, further, a 
high quality sinusoid, suggesting that the inferred rota- 
tion is not simply due to stochastic deviations from zero 
at a few position angles. The rms residual for the best- 



fitting model is 0.27, and the coefficient of determination 
R 2 = 0.97. 

Just 2.01% of our mock systems in set 1 had inferred 
rotation of matching or greater amplitude than that mea- 
sured for NGC 1846; however only about a tenth of these 
(0.23% of all systems) also had an rms residual value 
equal to or smaller than the best-fit model for NGC 
1846. A slightly larger fraction (0.39% of all systems) had 
matching or greater rotation amplitude together with an 
R 2 value equal to or greater than our best-fit model. 
Considering the lower bound on our NGC 1846 rotation 
amplitude of ~ 0.7 kms^ 1 , defined by the uncertaintly 
on the measurement, the fraction of mock systems with 
matching or greater amplitude rose to 8.88%. Again, 
however, only a small subset of these also had an equiv- 
alently high-quality rms value (1.31% of all systems) or 
R 2 value (0.97% of all systems). 

To give a visual indication of the goodness-of-fit crite- 
ria, we show in Figure [TT] several examples (taken from 
set 17, see below) where a rotation amplitude greater 
than that measured for NGC 1846 is inferred, but where 
the rms and/or the R 2 values indicate a poor quality si- 
nusoid. It is quite clear that these systems would not be 
mistaken for one in which a strong rotation signal had 
been reliably detected. We also show one example where 
stochastic fluctuations have resulted in a rotation curve 
that passess all our tests, and would have led to a false 
positive detection of rotation in the cluster. Fortunately, 
as outlined in Table 01 this is a rare occurrence. 

For each random realization of our sample we also cal- 
culated the mean velocity dispersion using the maximum 
likelihood technique described previously. We then found 
the average dispersion across all random realizations in a 
given set, along with the standard deviation in this value. 
The mean systemic velocity for each random realization 
is a natural by-product of the maximum likelihood cal- 
culation, and for completeness we determined the aver- 
age value for this quantity across each set as well. The 
results of this process are visible in Table 21 For the con- 
trol set 1, we recover an average systemic velocity and 
velocity dispersion precisely matching those values de- 
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rived for NGC 1846, including the uncertainties. This is 
to be expected, since the control run is the trivial case of 
re-measuring systems generated according to the velocity 
dispersion profile seen in Figure fTOl with no contribution 
from binary stars. 

We now consider the results for sets 2—17, which pos- 
sess binary population characteristics as outlined in Ta- 
ble HJ As expected, the additional velocity components 
due to the orbital motions of binary pairs do have some 
influence on our kinematic measurements. However, the 
effects do not appear to be strongly dependent on the 
adopted distributions for the orbital period or eccentric- 
ity. Changing the distribution for the mass ratio, q, has 
a more significant effect; however this is most likely due 
to the resulting change in the inferred binary fraction in 
the cluster - more binaries imply a stronger influence on 
the kinematic measurements. 

In terms of the inferred rotation, an increased percent- 
age of the mock systems have amplitudes of matching or 
greater amplitude than our NGC 1846 measurement. For 
fb = 0.33 these values lie in the range ~ 5.7—6.3%, while 
for h = 0.55 they are higher, ~ 7.7 — 8.3%. As before, 
however, only a very small subset of these systems pos- 
sess high quality sinusoidal rotation curves - the overall 
fractions are less than 1% for both goodness-of-fit crite- 
ria in all cases. Once again considering the conservative 
case defined by A lot ~ 0.7 kms -1 , the percentages all 
rise; however the overall fraction of systems in each set 
with a high amplitude and high-quality rotation curve is 
again always close to or below 1.5%. 

The mean kinematics of each set of models are inter- 
esting. Binaries have no effect on the average systemic 
velocity, which always comes out at the measured NGC 
1846 value. This is because the extra binary compo- 
nent of an individual radial velocity is equally likely to 
be pointing towards, or away from, the observer. How- 
ever the binary stars do appreciably inflate the mean 
cluster velocity dispersion. In the case where fb = 0.33 
the mean dispersion goes from 1.8 kms -1 to nearly 2.2 
kms , while for /& = 0.55 the dispersion is typically 
just above 2.3 kms . As with the rotation measures, 
these values are not strongly dependent on the adopted 
distributions for the orbital period or eccentricity. 

In summary, we draw the following conclusions from 
our Monte Carlo simulations. First, while it is not im- 
possible for rotation signals comparable to that which 
we measured for NGC 1846 to arise stochastically, it is 
a very unusual occurrence with < 1% of systems in any 
given set exhibiting such characteristics. The presence 
of unresolved binary stars increases the number of sys- 
tems for which a high amplitude of rotation could be 
inferred; however the rotation curves for such systems 
would invariably not be mistaken for the expected sinu- 
soidal variation of velocity with position angle. Thus, 
while we cannot rule out that the rotation we have ob- 
served in NGC 1846 is simply down to statistical fluke, 
we are better than 99% confident that it is a genuine 
signal. 

Second, it is likely that our measured velocity disper- 
sion for NGC 1846, a c i, along with our inferred central 
velocity dispersion ao, are both inflated above their true 
values due to the presence of unresolved binary stars 
in our sample. In our simulations, an input value of 
<7 c f = 1.81 kms -1 is increased to at most sa 2.33 kms . 
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Figure 11. Example rotation curves for three mock systems in 
set 17. We selected these to have comparable amplitudes to that 
measured for NGC 1846 (indicated, as previously, by the shaded 
region). The upper two panels show systems for which the best-fit 
models are of lower quality than than for NGC 1846, as defined by 
the rms of the fit and the coefficient of determination R 2 (see text). 
In the lower panel we show one rare example of a high quality curve 
arising purely by chance. This particular realization had the best 
quality fit of all 10 J mock systems in set 17. 

Since the various sources of dispersion add in quadrature, 
this implies the true a c i could be as low as~ 1.10 kms -1 , 
and the true ao as low as s» 2.05 kms -1 . By comparison, 
our typical measurement uncertainty of « 0.6 kms -1 
on each individual velocity inflates the true dispersions 
by less than 0.1 kms -1 (note that this contribution is 
accounted for by the maximum likelihood technique we 
employed to measure the cluster dispersion). 

5. DISCUSSION 

The key point of interest arising from our analysis of 
the internal kinematics of NGC 1846 concerns the de- 
gree of systemic rotation, which appears high compared 
with, for example, many Milky Way globular clusters, 
or indeed the few Magellanic Cloud clusters for which 
such measurements exist. We find a ratio A rot /ao = 
0.44 ± 0.16 if we adopt the mean amplitude of rotation 
A ro t = 1.1 ± 0.4 kms -1 . However, we also demonstrated 
that the maximum amplitude of rotation in NGC 1846 
might well be as high as ps 2 kms -1 (especially since 
there is the additional unknown factor sinz required to 
deproject our rotation measurement - note that the mean 
value of sin i assuming a uniform distribution of inclina- 
tion angles is 2/n). In this case, A rot /cro ~ 0.8. We 
also showed that unresolved binary stars in our sample 
probably inflate the measured velocity dispersion by up 
to pa 0.5 kms -1 , implying that the fraction of ordered 
motion with respect to pressure support is likely to be 
even larger still, A Iot /ao <; 1-0. 

These latter two estimates are commensurate with the 
highest ratios measured in the sample of 24 Galactic 
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globular clusters compiled by iBellazzini et al.l (|2012f ) ; in- 
deed even our lower limit of A rot /cro = 0.44 would place 
NGC 1846 in the top 25% of this sample. They are 
also somewhat larger than t he values measured for young 
Magellanic Cloud clusters. iFischer et al.l (|1992a[) found 
A mt /a = 0.45 ± 0.20 in the sa 1 00 Myr old cluster NGC 
1866, while iFTscher et all ()1993f ) measured a comparable 
degree of rotat ion in the ~ 50 Myr old cluste r NGC 1850. 
More recently, Hcn ault-Brunet et al.l (2012a. b) measured 
-Arot/fo = 0.60 ±0.30 for the very young (sa 3 Myr) mas- 
sive cluster R136 at the centre of 30 Doradus. Notably, 
however, the only other intermediate-age LMC cluster in- 
vestigated on a star-by-star basi s, NGC 1978, shows no 
significant evidence for rotation ()Fischer et al.lll99 2b). 

When comparing the degree of internal rotation in 
NGC 1846 with rotation seen in Galactic globular clus- 
ters, it is important to bear in mind that NGC 1846 
is still dynamically quite young. As demonstrated in 
Section 14.41 its central and median relaxation times are 
greater than the cluster age. It is possible, even likely, 
that these time-scales were much shorter early on in th e 
cluster's evolution (see e.g.. IMackev et al.l [20071 12008bf) : 
however it is clear that in the ancient Galactic globular 
clusters any ordered motions have had far more oppor- 
tunity to become randomized, thus surpressing A rot /<7o 
relative to NGC 1846. 

This is not the case for the young LMC clusters NGC 
1850 and 1866, and R136, which, li ke NGC 1846 should 
be dynamically unevolved (e.g., IFischer et al.l Il992al 
119931: iHenault-Brunet et alJl2012a[ ). That the degree of 
rotation observed in NGC 1846 is, arguably, greater than 
in these young systems is a striking result. If, as we sus- 
pect, A rot /ao « 0.8 — 1.0 then ordered rotation is com- 
parable in importance to random motions in providing 
support a gainst gravitational collapse. Fo llowing the dis- 
cussion in Hcnault- Brunet et al.l (|2012al ). for this range 
in A Iot /ao, between sa 35 — 45% of the total cluster ki- 
netic energy would be in rotation. 

NG C 1846 is only mildly elliptical - iGoudfrooii et al.l 
(2009) measure e = 0.12 ± 0.02. However, a cluster's el- 
lipticity is not neces sarily a good indicatio n of its degree 
of internal rotation. IBel lazzini et"aT1 (|2012h found no cor- 
relation between these two quantities in their sample of 
Galactic globular clusters. Furthermore, the lone globu- 
lar cluster associated with the isolated Local group dwarf 
irregular galaxy WLM, shows no evidence for strong in- 
ternal ro tation despite display ing a high degree of el- 
lipticitv ([Stephens et alJ l200ffl. NGC 1978, for which 
IFischer et alj (|1992bl ) found no evident signs of rotation, 
al so has a very fla t tened shape with e w 0.3. 

Bcllazzi ni et al.l (|2012ft demonstrated that for Galac- 
tic globular clusters A lot /ao correlates quite strongly 
with [Fe/H], such that more metal-rich systems have 
stronger internal rotation. With [Fc/H] sa —0.4, NGC 
1846 fits this correlation well, although we again cau- 
tion that we are comparing dynamically evolved systems 
wit h a dynamically young system. The suggestion made 
bvlBeTI azzini et al.l ([20121 ) is that the observed correla- 
tion between [Fe/H] and A rot /ao may hint that dissipa- 
tive gas dynamics plays a significant role in the process 
of cluster formation, due to the fact that a larger metal 
content in a gas should imply a higher efficiency in energy 
dissipation via atomic transitions. 

This idea is particularly relevant to the present case, 



where we strongly suspect t he pres e nce of multi ple stellar 
generations in NGC 1846. iBekkil (pOlOl |201~1 uses hy- 
drodynamic simulations to investigate the formation of a 
second generation of stars at the centre of a globular clus- 
ter following the accretion of ga s (including AGB ejecta ) 
into its potential well (see also iBekki fc Mack cv 2009) . 
He finds that if the first generation of stars possesses even 
a very small net angular momentum, gaseous dissipation 
during accretion onto the cluster center leads to a dy- 
namically cold, rapidly rotating second generation with 
j4rot/co > 0.8. The second generation is initially very 
centrally concentrated but this nested structure, and the 
strong rotational signature, ought to be smoothed out 
as a result of relaxation processes during the subsequent 
long term evolution of the cluster. Since NGC 1846 is, 
dynamically speaking, still quite young, both the central 
concentration of the yo unger generation (as measured by 
iGoudfrooij et al.ll2009j ). and the rotational kinematics (as 
seen in the present work), are apparently still evident. 
Our possible detection of increasing rotation towards the 
cluster center is also consistent with this picture. 

A more sensitive test of the model described above 
would be achieved if we could split our kinematic sample 
into earlier and later generations; however, to presage the 
results of our chemical abundance analysis somewhat (see 
Paper II), there is no obvious marker for achieving this 
- and in any case the present ensemble is probably too 
small for useful subdivision. The most reliable method 
would be to target a large number of stars across the 
EMSTO in order to directly correlate kinematics against 
position on the CMD and within the cluster. Such mea- 
surements, if at all possible, would require a considerable 
investment of telescope time but would provide critical 
information to test the different formation hypotheses. 

Nonetheless, the fact that we have detected strong 
rotation in the first EMSTO cluster subjected to de- 
tailed dynamical examination is suggestive that this 
may be an important feature of these systems. In fu- 
ture it will be critically important to enlarge the sam- 
ple of intermediate-age Magellanic Cloud clusters for 
which internal kinematics have been measured, includ- 
ing both systems with and without an EMSTO. Real- 
istic TV-body modelling of clusters such as those pre- 
dicted by hydrodynamic simulations, in which there is 
a centrally-concentrated, strongly-rotating second gen- 
eration of stars embedded in a more diffuse, pressure 
supported first generation, would also be extremely use- 
ful (although we recognize the current limitations placed 
on such models due to the maximum particle number of 
~ l-2x 10 5 ). An issue of particular interest is how the 
spatial and kinematic distinctions between the different 
generations propagate through the dynamical evolution 
of the system, especially where the cluster is strongly 
mass segregated at early times such that violent relax- 
ation due to rapid stellar mass-loss drives significant ex- 
pansion, as appears to be n ecessary for the in termcdiate- 
age EMSTO clusters (e.g. . Keller et^[20lll ). 

6. SUMMARY & CONCLUSIONS 

In this paper we have described a detailed set of 
VLT/FLAMES observations of red giant stars in the 
peculiar intermediate-age LMC star cluster NGC 1846, 
along with the data reduction procedure we employed 
to extract and process individual spectra. In total, we 
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targeted 29 stars within the nominal boundary of NGC 
1846, of which 21 possess radial velocities indicating their 
membership of the cluster at high confidence. In addi- 
tion, we targeted the planetary nebula Mo-17, and the 
radial velocity of this object indicates that it too is a 
member of the cluster. 

We have used our spectra to investigate the elemen- 
tal abundance patterns present in NGC 1846, including 
the possibility of star-to-star variations in light element 
abundances. These results will be presented in a forth- 
coming work (Paper II). In the present paper we took 
the radial velocity measurements for our sample and used 
these to conduct a thorough analysis of the internal kine- 
matics of the cluster, with the following results: 

• NGC 1846 exhibits a significant degree of systemic ro- 
tation. The mean amplitude is A Iot = 1.1 ±0.4 kms -1 , 
with the rotation axis oriented at 60±20° east of north. 
There are indications that the rotation signal may vary 
with position in the cluster such that the amplitude in- 
creases towards the center and peaks somewhere within 
the half-light radius. The maximum amplitude may 
well be as high as A Iot s» 2 kms -1 , especially consid- 
ering that there is also a correction factor sin i for the 
unknown inclination of the rotation axis to the line 
of sight. An extensive suite of Monte Carlo models 
suggests that, because of the relatively small size of 
our sample and the presence of a significant popula- 
tion of unresolved binary stars in NGC 1846, stochas- 
tic fluctuations could reproduce the observed rotation 
curve; however this only occurs very rarely - less than 
« 0.3% of the time if no binaries are present, or less 
than sa 1% of the time for a cluster binary fraction of 
up to fb — 0.55. 

• We measure a mean velocity dispersion a c i = 1. 81^^29 
kms -1 . Assuming a simple parametrization of the ve- 
locity dispersion fall-off with radius, the implied cen- 
tral velocity dispersion in the cluster is en. = 2.52+°'^ 
kms -1 . Our Monte Carlo modeling suggests that the 
presence of unresolved binary stars in our sample could 
substantially inflate these quantities. If the binary 
fraction /& = 0.55, the true values could be as low 
as a c i ~ 1.1 kms -1 and ctq ~ 2.0 kms -1 . 

• The ratio of ordered motion to pressure support is for- 
mally A rot /(io = 0.44 ± 0.16; however, accounting for 
the probable maximum amplitude of rotation in the 
cluster, the inclination factor sini, and the contribu- 
tion of binary stars to inflating the observed veloc- 
ity dispersion, this quantity is likely to be as high as 
0.8 — 1.0. In this case, between « 35 — 45% of the total 
cluster kinetic energy would be in rotation. 

• Under the assumption that mass follows light in the 
cluster, the mass of NGC 1846 is (8.4±J-j) x 10 4 M ? and 
the implied mass-to-light ratio is 0.5915'Jq, consistent 
with predictions made purely on consideration of its 
consitutent stellar populations. Note, however, that 
these quantities are derived assuming a "dispersion- 
only" cluster (Eq. HJ. If internal rotation provides an 
important contribution against gravitational collapse, 
as seems probable, more sophisticated modelling will 
be required to obtain reliable estimates of mass and 
M/L v . 



• The median relaxation time for NGC 1846 is t r h = 
2.61^2 Gyr, indicating that the cluster is dynamically 
youthful. Hence any kinematic signatures encoded 
during its formation ought to remain present. 

The observation that substantial rotation is present in 
NGC 1846, at a magnitude comparable to that of the 
velocity dispersion, is consistent with the predictions of 
simulations modeling the formati on of multip l e gen era- 
tions in globular clusters (see e.g., lBekOl2lM[20llh . It 
would be of significant interest to improve our knowledge 
of the internal kinematics of this cluster by extending 
the present work to a much larger sample, ideally one 
in which the multiple generations could be easily identi- 
fied. Similarly, by extending our analysis to additional 
intermediate-age Magellanic Cloud clusters, both with 
and without the extended main-sequence turn-off mor- 
phology, we could hope to learn whether strong internal 
rotation is a key signature of the formation of clusters 
with multiple constituent stellar populations. 
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